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Zero-temperature or quantum phase transitions in itinerant electronic systems both with and 
without quenched disordered are discussed. Phase transitions considered include, the ferromag- 
netic transition, the antiferromagnetic transition, the superconductor-metal transition, and various 
metal-insulator transitions. Emphasis is placed on how to determine the universal properties that 
characterize these quantum phase transitions. For the first three of the phase transitions listed 
above, one of the main physical ideas established is that in zero-temperature systems there are 
soft or slow modes that exist in addition to the soft order parameter fluctuations, and that these 
modes can couple to the critical modes. These extra soft modes are shown to have a profound ef- 
fect on the quantum critical properties. For quantum phase transitions involving zero wavenumber 
order parameters, i.e., the ferromagnetic and superconductor-metal transitions, these extra modes 
effectively lead to long-ranged effective interactions between order parameter fluctuations, which in 
turn lead to exactly soluble critical behaviors. For the antiferromagnetic case, we argue that while 
in low enough dimensions disorder fluctuation effects tend to destroy long-range order, quantum 
fluctuations counteract this effect and in some parameter regions manage to re-establish antiferro- 
magnetic long-range order. For the metal-insulator transition, some recent new ideas are reviewed. 
In particular, it is pointed out that for interacting disordered electrons, one expects that in high 
dimensions the metal-insulator transition is related to the phase transition that occurs in random- 
field magnets in high dimensions. If the analogy also holds in three dimensions this suggests that 
the metal-insulator transition might have glassy characteristics. 
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I. INTRODUCTION 



The study of continuous phase transitions EJ has led to great advances in condensed matter physics. The insights 
gained from these studies ||l| have had remarkably broad effects, which is somewhat suprising given that phase 
transitions occur only at very special points in the phase diagram. In particular, phase transition theory spawned 
the application of renormalization techniques in condensed matter physics |^ , which have proved to be very powerful 
in a more general context. Wilson's breakthrough paved the way for the understanding of thermal or classical phase 
transitions, which occur at non-zero temperature and are driven by thermal fluctuations. Let the critical temperature 
be Tc > 0, and denote the distance from the critical point hy t = l — T/Tc- Then it turns out that the critical behavior 
asymptotically close to the critical point is entirely determined by classical physics. This can be seen as follows. Upon 
approaching the critical point, there is a diverging characteristic length in the system, viz. the correlation length 
f ~ 1^1 with 1/ an exponent that is characteristic of the class of phase transitions under consideration. Also, the 
characteristic time for relaxation towards equilibrium diverges like ~ ^'^j with z another exponent. This slow 
relaxation is caused by slow fluctuations of some thermodynamic quantity, the order parameter, that is characteristic 
of the transition. These fluctuations relax on microscopic time scales far from the transition, but become infinitely 
slow as t — *■ 0. This means there is a characteristic frequency, ujc ~ that vanishes at criticality. Now quantum 

effects will be unimportant as long as huic << ksT, which is always the case sufficiently close to a critical point with 
a non-zero Tc. As a result, the critical behavior at, say, the Curie point in Iron, the A-transition in liquid Helium, or 
the superconducting transition in Mercury, are determined entirely by classical physics, notwithstanding the fact that 
quantum mechanics is crucial for the microscopic mechanism underlying the transition and, in the last two examples, 
for the properties of the ordered phase. 

These considerations raise the interesting question of what happens when we follow a phase separation line in a 
phase diagram down to zero temperature. For definiteness, let us consider a ferromagnet, or a superconductor, where 
we decrease by diluting it with a non-magnetic or non-superconducting material. Then we obtain a schematic 
phase diagram like the one shown in Fig. |l|, with J some parameter that triggers the transition. Normally, J is fixed 
for a given material, and one observes the transition by lowering the temperature. However, if we imagine crossing 
the phase separation line at T = by varying J, then the above arguments suggest that quantum effects should 
be important for the transition, since Tc = 0. We therefore expect a critical behavior in this case that is different 
from the one observed when crossing the phase separation line at any non-zero T. This question is not academic, 
as was first pointed out by Suzuki |^ and Hertz If one triggers the transition at low but non-zero T, then 
one observes a 'crossover' from a region that is characterized by quantum critical behavior to one showing classical 
critical behavior, as indicated in the figure. Moreover, there are phase transitions that occur only at T = 0, the most 
important examples being various metal-to-insulator transitions, which are studied experimentally by observing a 
smeared transition at very low temperatures. 

Such phase transitions that occur in quantum mechanical systems at T = as a function of some non-thermal 
control parameter are called quantum phase transitions (QPTs). Like their finite-temperature counterparts, which 
for distinction are often referred to as thermal or classical phase transitions, they are characterized by a diverging 
correlation length ^, and a diverging relaxation time ^r- However, the critical fluctuations that lead to these diverging 
length and time scales are quantum fluctuations rather than thermal ones. Apart from the origin of the relevant 
fluctuations being different, quantum phase transitions also differ from their classical counterparts in that their 
dynamic and static critical behaviors are coupled together. Technically this implies that one cannot solve for the 
thermodynamic properties of the phase transition without also solving for the dynamics. 

There are a number of ways to see the coupling between statics and dynamics at QPTs. From a general scaling 
viewpoint one argument is as follows. It is obvious from Fig. |l| and the accompanying discussion that increasing the 
temperature from zero takes one away from the quantum critical point, much like increasing the frequency from zero 
takes one away from either a quantum or a classical one. Now the fact that Planck's constant h and Boltzmann's 
constant fcs are choajjn arbitrarily implies that in fundamental units, frequency, energy, and temperature all have the 
same dimensions. □ From a scaling viewpoint, this implies, for example, that if a system is at a quantum critical 
point, either going to a finite frequency or increasing the temperature will affect the system in the same way, and be 
characterized by the dynamical scaling exponent that describes the divergence of £,r and is usually denoted by z. Since 
the thermodynamic observables certainly depend on the temperature, this further implies that the scaling relations 
or homogeneity laws for these quantities will also depend on z. That is, the static and dynamic critical behaviors are 



^We will deal only with continuous transitions, and will simply refer to them as 'phase transitions' or 'critical points'. 
■^This has recently been emphasized in Ref. pi. 
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FIG. 1. Schematic phase diagram with J some parameter that is responsible for the ordering. The quantum critical point 
is at (T = 0, Jc). The dashed lines denote the boundaries of the critical region, and the dotted lines mark the crossover from 
quantum to classical critical behavior. The regions labeled by c and q show classical and quantum critical behavior, respectively. 



coupled together. A more mathematical argument is that the standard Matsubara frequency technique implies that 
frequency and temperature are intimately related and should scale in the same way. 

From a more technical point of view, a classical partition function involves separate integrals over momenta and 
positions that are decoupled since the kinetic energy depends only on the momenta, and the potential energy only on 
the positions. One can therefore solve for the thermodynamics without having to deal with the dynamics explicitly. 
In cuiantum statistical mechanics, on the other hand, the partition fimction of a system with Hamiltonian H is given 

as a 

Z = Tre^-^/'^'^^ . (1.1) 

All methods for calculating the above trace make use of the imaginary time technique |^ in one form or another. 
The latter is based on the observation that e~^^ , with p = l/fc^T, is formally identical with the time evolution 
operator e'^*/'', with the time t replaced by the imaginary quantity ihp. One then proceeds to divide the imaginary 
time interval [0, /3] into many infinitesimal subintervals for which the time evolution can be calculated. The coupling 
between statics and dynamics then arises from the fact that the Hamiltonian taken at some imaginary time, H{t), 
does not commute with H taken at another imaginary time r'. 

Although QPTs are currently a topic of great interest in condensed matter physics, they are actually a subject with 
a long history. The earliest work all dealt with quantum spin systems, and especially one-dimensional chain models. 
Important work was done by Bethe, Bloch, Tomonaga, and Lieb. In the mid-1970's the QPT problem was cast into 
a more general theoretical framework. Suzuki |^ used the so-called Trotter formula, which is one way to rewrite 
the trace in Eq. (^^) by dividing [0, (3] into sub-intervals, to map some d-dimcnsional quantum spin problems onto 
(d + l)-dimensional classical statistical mechanical problems. The extra dimension was for imaginary time, which is 
of infinite extent in the zero-temperature limit (/3 — > oo), effectively increasing the dimension of the system. Beal- 
Monod noted that this mapping made it obvious how to generalize the renormalization group techniques that were 
being used to describe classical phase transitions for an application to quantum transitions. Motivated by these works, 
Hertz [Q made a number of important advances. First he showed how to derive Landau-Ginzburg- Wilson (LGW) 
order-parameter functionals, starting from a fermionic description of itinerant electron systems that writes the trace 



^We will denote by Tr traces over all degrees of freedom, including continuous ones, and by tr traces over discrete degrees of 
freedom that are not shown explicitly. 
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in Eq. (1.1) as a functional integral. As expected, he showed that to describe QPTs one needs to consistently describe 
both space and (imaginary) time fluctuations together. He further showed that Suzuki's mapping of a d-dimensional 
quantum system on a (d + l)-dimensional classical one was too restrictive for generic systems. He found that, more 
generally, a d-dimensional quantum system is related to a (c? + z)-dimensional classical one, with z the dynamical 
scaling exponent for the quantum phase transition. The need for z ^ \ arises since in general in a condensed matter 
system, space and time do not scale in the same way. For instance, in a system with diffusive dynamics, one expects 
frequency to scale like the wavenumber squared, which suggests z = 2. The special case z = 1 occurs when space and 
time happen to scale in the same way. Hertz then showed that the same renormalization group ideas that had been 
used to characterize classical phase transitions could be applied to the quantum case for arbitrary values of z. 

Hertz concluded that most QPTs in three-dimensional (S-d) systems are trivial from a phase transition viewpoint 
in the sense that they are characterized by Landau or mean-field critical exponents. This conclusion was based on 
the observation that the mapping d —> d + z, with z > 1, lowers the upper critical dimension djj", above which the 
transition is correctly described by a mean-field/Gaussian theory, from djj" = 4 for most classical phase transitions to 
d+ = 4 — z < 3 for the quantum transition. For systems with z = 1 the critical behavior would still be mean-field 
like up to logarithmic corrections to scaling. The only exceptions would be low-dimensional (d < 2) systems. 

One of the main points of the current review is to emphasize that most QPTs are more complicated than suggested 
by the above arguments. These complications arise because in many quantum systems there are soft or slow modes that 
exist in addition to the critical order parameter fluctuations. These "extra" modes arise either from conservation laws 
(recall that statics and dynamics are coupled at QPTs), or from Goldstone modes that result from some spontaneously 
broken continuous symmetry. These modes will in general couple to the critical modes and influence the critical 
behavior. If one insists on integrating out these degrees of freedom in order to obtain a standard LGW theory that 
is formulated in terms of the order parameter only, then this LGW theory becomes non-local. Physically, the extra 
modes lead to long-range interactions between the order parameter fluctuations, which in turn leads to a non-local 
LGW functional. This breakdown of Hertz's theory becomes particularly obvious in systems with quenched disorder. 
In the cases considered by Hertz, he again found that the correlation length exponent had its mean-field value, 
V = 1/2, for all physical values of the dimensionality. However, it is now known (though it was not known when Hertz 
derived his theory) that v must obey the inequality v > 2/ din generic disordered systems This implies that the 
mean-field result cannot be correct for d < 4. 

The LGW theories mentioned above are one example of what is called effective field theories for condensed matter 
systems. The main motivation for deriving and studying such effective field theories is that often one is not interested 
in all of the detailed behavior that arises from the microscopic properties of a system. For instance, one might 
ask about universal effects that are the same for large classes of systems, even though the members of these classes 
are different with respect to their composition, their band structure, and other microscopic details. Such universal 
properties arise from long-wavelength and low-frequency fluctuations, i.e. processes that occur on length and time 
scales that are large compared to the microscopic lengths and times that determine, say, the band structure. This 
in turn implies that only the dominant fluctuations in this limit need to be retained explicitly and treated carefully 
in order to describe the universal properties. Of particular interest are soft modes, i.e. excitations that are massless 
in the limit of zero frequencies and wavenumbers. All other degrees of freedom can be integrated out, often in 
a crude approximation, without affecting the description of the universal behavior. The resulting theory for the 
most important fluctuations in a specified system is called an effective field theory. This is quite analogous to the 
situation in high-energy physics, where one also deals with effective field theories that describe the low-energy or 
long-distance behavior of some underlying microscopic theory. The main difference is that in condensed matter 
physics the microscopic theory is known, but too complicated to solve in general, and the problem consists of deriving 
an effective theory that is simple enough to extract the behavior of interest, while in high-energy physics the effective 
theory is known (i.e., it has been guessed) and the problem consists of inferring from it a more microscopic theory. 
In high-energy physics one thus proceeds from low to high energies, while the opposite is true in condensed matter 
physics. The most well-known examples of effective field theories in a condensed matter context are those that are 
used to describe the universal behavior near critical points. The concept is much more general, however. Even far 
from any phase transitions, soft modes will arise from either the conservation laws, or from the existence of Goldstone 
modes in a broken symmetry phase. Both of these mechanisms are common in quantum statistical mechanics, and 
the resulting soft modes determine the physics at large distances and low energies. In addition to these soft modes, 
near a phase transition the critical order parameter fluctuations are also soft. One of our objectives is to review how 
to derive effective fleld theories, both near and far from phase transition lines or points. 

The current chapter is organized as follows. In Section^ we give a brief discussion of scaling theory near (continuous) 
quantum phase transitions. We start Section [II with a general discussion of the functional field theoretic description 
of itinerant electronic systems. We then discuss how to derive effective field theories for the relevant degrees of 
freedom, and show how to describe soft modes in general disordered electronic systems. We begin Section IV by 
reviewing the ferromagnetic phase transition in both disordered and clean itinerant electron systems. We will see that 
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weak localization effects, and their clean system analogs, determine the universal properties of this phase transition. 
We then examine the antiferromagnetic transition in disordered itinerant systems. The most important result is that 
quantum fluctuation effects stabilize antiferromagnetic long-range order in certain parameter regions. In Section ^ 
we discuss some recent work on the superconductor-metal transition in bulk [d = 3) disordered electron systems. In 
Section^ we consider the disordered Fermi liquid phase, and its instability against metal-insulator transitions of the 
Anderson-Mott type. First we briefly review the theory based on Finkel'stein's generalized nonlinear sigma-model, 
and then we discuss a recent order-parameter approach to the Anderson-Mott transition. We discuss similarities 
to the classical random-fleld magnet problem and discuss possible implications for the Anderson-Mott transition in 
three dimensions. 

There are other interesting and important examples of quantum phase transitions that we will not discuss, for 
instance quantum spin chains, quantum spin glasses, and phase transitions in quantum Hall systems. A brief review 
of the latter subject has recently been given by Sondhi et al. §. 



II. SCALING AT QUANTUM CRITICAL POINTS 

Scaling ideas have been used with great success in the description of continuous thermal phase transitions since 
the mid-1960's [Q. These same ideas, properly reinterpreted, can be applied to the QPT problem. Since the scaling 
description of quantum phase transition has been recently reviewed in a number of places P,p^ we will be brief here. 

In general, the mapping d — s- d + z, discussed in Sec. ^ above, suggests that the classical scaling relations are valid 
at the QPT, with the only change being the above substitution. Although this is the standard lore, and it is true in a 
formal sense, we will repeatedly encounter situations where it is not valid. When such a breakdown of general scaling 
happens, we will mention it explicitly. More specifically, our main conclusions will be as follows. (1) The classical 
scaling relations that are independent of dimensionality are also valid at QPTs. (2) The classical scaling relations that 
do depend on d ('hyperscaling relations') are formally valid with the substitution of d — > d + z. However, since many 
of the important QPTs are above their upper critical dimensionalities, these scaling relations in general break down 
because of the presence of dangerous irrelevant variables. □ Further, the situation is much more complicated than 
in classical systems above their upper critical dimensionality because, as mentioned in Sec. ||, the LGW functional 
for describing QPT's is often a non-local one. This feature has been shown to lead to nontrivial, dimensionality 
dependent, exponents and scaling relations, even above the upper critical dimension for the QPT. (3) There are 
scaling relationships at QPTs that have no classical analog. These identities follow from the scaling equivalence of 
energy, inverse time, and temperature. They depend on the dimensionality and can also be invalidated by dangerous 
irrelevant variables. 

All scaling theories for continuous phase transitions start with the idea of a divergent correlation length scale, ^. If t 
is the dimensionless distance from the critical point, then the divergence of ^ is characterized by the critical exponent 
I', 

i - tV • (2-la) 

Typically t is defined to be greater than zero in the disordered phase, and less than zero in the ordered one. At 
thermal phase transitions, t is given by the reduced temperature, and it is very easy to control experimentally. For 
QPTs t might represent the distance from a critical interaction strength, or from a critical impurity concentration. In 
the quantum case, t is therefore not so easily controlled, and to vary t different samples are usually used. □ In general 
this means that experimentally probing the quantum critical region, with the precision necessary for measuring critical 
exponents, is very difficult. 

As we have mentioned in Sec. |, in addition to the diverging correlation length we need to deal with a diverging 
time scale. As the critical point is approached from, for example, the disordered side, the order that is being built up 
decays on a longer and longer time scale. This correlation time i^^- diverges at the transition, and the singularity is 
characterized by the dynamical scaling exponent z, 

C.^r • (2.1b) 



*The concept of dangerous irrelevant variables is defined at the end of Sec. ||. 

^An exception is the stress tuning technique, where uniaxial or hydrostatic stress is applied to a sample that is very close to 
the transition. The resulting change in the microscopic structure drives the system across the phase separation line. See Ref. 
[hll for an application. 
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Equations ( |2.1| ) represent an assumption, namely that both ^ and diverge as some power of t. There is no principal 
reason why the singularity might not be stronger or weaker than that, and we will indeed encounter examples of 
exponential divergencies in Sees. ^ and VI below. For now we will continue, however, with the discussion of the more 



conventional power-law scaling assumption. 

Before formalizing the scaling ideas, we will first give two examples of how scaling (or hyperscaling) is used. First 
consider the singular part of the free energy density, fs . The dimension of fg is that of an inverse volume times an 
energy. Assuming that f and are the only relevant length and energy □ scales in the problem, this suggests 

/s^r''+"^ . (2.2a) 



Note that for the classical problem, Eq. ( p. 2a ) is obtained with d + z replaced by d , since in that case the relevant 



energy scale as ^ ^ oo is just the critical temperature, Tc ^ 0. Using Eq. ( 2.1b ) with ^ ^/T yields for the 
temperature dependence oi fs &t t = 

/,(t = 0,r)-Ti+^/^ . (2.2b) 

The singular part of the specific heat thus behaves as 

C - T"^/' . (2.3) 

That is, by using scaling ideas one can express the critical behavior of C in terms of the dynamical critical exponent. 

As another example, consider the conductivity, a. In fundamental units, setting h — = = 1 (with e the 
electron charge), the dimension of a is an inverse length to the power (d—2). At a phase transition where a vanishes 
in the 'disordered' phase and is nonzero in the 'ordered' phase, i.e., at a metal-insulator transition, u this suggests 
that a vanishes as 

a{T = 0, t) - ^-(''-2) _ ^2.4a) 
We see that the conductivity exponent s, defined by cr ^ i*, is related to the correlation length exponent v by 

s = iy{d-2) . (2.4b) 



Equation (2.4b) is called Wegner's scaling relation |12| ]. 

The above scaling arguments have used only what in the literature is called a naive or engineering (with apologies 
to engineers) dimensional analysis. For instance, we have assumed that since the conductivity is measured in units of 
centimeters to the power 2 — d, it will scale like ^2^**. In general this type of argument is not correct. The conductivity, 
for example, could acquire an 'anomalous dimension' 9 by scaling like a ^ ^s-d+e ^-e^ where the 'wrong' dependence 
on the length scale ^ is compensated for by a dependence on some microscopic length scale a. However, for some 
quantitites, such as / and a above, there are additional arguments that show that they have zero anoma lous d imension 



so that the above argument is correct. There is, however, yet another way for scaling relations like Eq. (2.4b) to break 
down, which we will discuss at the end of the current Section. 

Let us now formalize the scaling approach a bit. Suppose we are interested in an observable Q as a function of t, T, 
frequency Q, wavenumber fc, and possibly an external field h. The power law singularities that characterize continuous 
phase transitions imply that the singular or scaling part Qs of the observable is a generalized homogeneous function of 
the appropriate variables. That is, the ususal scaling ansatz, which typically can be derived by using renormalization 
group methods, is 

Qs{t,T,h;k,n) = b-^'^^Qs{b^^''t,¥T,by'^h;bk,b''n) . (2.5) 
Here b is an arbitrary length resc alin g factor, and the number [Q] is called the scale dimension of the quantity Q. 



The remaining exponents in Eq. (2.5) have been chosen so that b really is the length rescaling factor. For example, 
it Q — £, is the correlation length, then atT = /i = fc = r2 = 0, and using [£] = —1 (as one should if 6 is a length 
rescaling factor), we have ^{t) — b^{b^/'^t). Choosing b = l/t" gives ^{t) = ^,{1)/^, i.e., £, correctly diverges as t~'' as 
t — > 0, and ^(1) is the value of the correlation length far away from the critical point. 



^We will use units such that ft = 1 in what follows. 

''We will see in Sec. ^ in which sense the metallic phase can be considered the ordered one. 
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To illustrate the usefulness of these homogeneity laws, let us consider the conductivity a again, as a function of 
T, and 17, 



By choosing 6 = Eq. ( |2.6a| ) gives, 

and by choosing b — XjT^I^ we obtain, 

a{t, T, 17) = T^''-^)/^ a{tlT^l''\ 1, 17/r) 



(2.6a) 



(2.6b) 



(2.6c) 



Here we have defined the functions F„(x,y) = a(\^x^y) and Ga-{x,y) = a{x,l,y). The important point is that, say, 
at 17 = the quantity a/V'^'^^''^^ is a function of T /t'^'^ alone, rather then of t and T separately. The resulting collapse 
of data on a single scaling curve has been historically important for the confirmation of scaling at thermal phase 
transitions a nd it is equally useful at QPTs. We will consider an example in Sec. VI B 3 below. 

Equations ([2 . 6| ) reflect the fact that t is not the only relevant distance from the critical point. At zero frequency 
and wavenumber there are two relevant variables: t and T. For either t ^ Q ov T ^ Q the system is away from the 
critical surface, and the exponents v and z are measures of the relevance of t and T, respectively. Of course, in any 
system there are numerous other var iables that are irrelevant. Let us denote a generic irrelevant variable by u. In 



a homogeneity law such as Eqs. (2.6), these irrelevant variables also occur, but they have negative scale dimensions, 
[u] < 0. We generalize Eq. ( 2.6a ) at 17 = to 



Choosing 6=1 /f^ now gives 



cr(t, T, u) = b-'^'^-^^(7{b^^''t, b^T, fe^u) 



a{t, T, u) = r^'^-^) (7(1, ^/^^^ i-''["lu) 



which for f — > yields (remember [u] < 0) 

a{t 0, T, u) = t"^'^-^^ a(l, T/^''^ 0) 



(2.7a) 



(2.7b) 



(2.7c) 



We see that close to the critical point, t ^ 0, the dependence on the irrelevant variable drops out. This is a very 
important result and is one of the reasons for the universal behavior near continuous phase transitions. 

We conclude this subsection by introducing the concept of a dangerous irrelevant variable |Q. Let us use Eq. 
( 2.7c ) as an example. Implicit in the last equality in that equation is the assumption that a{l, T/f^^jU = 0) is neither 
zero nor infinity. This is a nontrivial assumption which can break down. If it does, then the irrelevant variable 
u is called a dangerous irrelevant variable and actually modifies the scaling behavior, and scaling equalities, even 
though it is irrelevant in the technical sense of the renormalization group. Whether or not a given irrelevant variable 
is dangerous with respect to a certain observable can usually not be determined on general grounds, but requires 
explicit calculations. In the context of thermal phase transitions, it is dangerous irrelevant variables that lead to a 
breakdown of hyperscaling above the upper critical dimensionality ||l| . In the following sections we will repeatedly 
see that dangerous irrelevant variables play an even more important role in Q PTs, m ostly due to the small value of 
t he up p er cr itical dimensionality of the latter. In particular, we will see in Sec. VI B 2 that Wegner's scaling law, Eqs. 



(2.4b), (2.7c), is violated at certain metal-insulator transitions due to a dangerous irrelevant variable. 



III. FERMIONIC FIELD THEORY 



The many-fermion problem has a long history, due to its importance with respect to electrons in condensed matter 
systems. In order to describe nonperturbative effects in these systems, such as phase transitions involving spontaneous 
symmetry breaking, it is convenient to go to a field theoretic description. In such a formalism the important underlying 
symmetries are most apparent, and, at least formally, it is relatively easy to integrate out the degrees of freedom that 
are irrelevant for describing the physics one wants to focus on. Here we will use standard field theoretic techniques 
to sketch the development of a theory for disordered interacting electrons, following along the lines of Ref. [p^. 
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A. Grassmannian field theory 



Since the description of fermions involves anticommuting variables, any field theory or functional Feynman path 
integral description for electrons must be formulated in terms of anticommuting or Grassmann variables Jl6t . For the 
sake of simplicity we consider a model for a homogeneous electron fluid subject to a random potential that models the 
quenched disorder. For most of the effects we want to study, the inclusion of microscopic details like the underlying 
crystal structure and the resulting band structure effects, or the explicit inclusion of phonons, is not necessary. The 
partition function can then be written Q 



Z = j Z?[V',V]exp[5] . (3.1) 
Here the functional integration is with respect to Grassmann valued fields, V' and ip, and the action S is given by 

S ^ ~ dx'^il)c{x)dri'o{x) + 5o + Sdis + Sint ■ (3.2a) 

•' a 

We use a {d + l)-vector notation, with x = (x, t), and J dx = Jy dx dr. x denotes position, t imaginary time, V 
is the system volume, /3 = l/T is the inverse temperature, and a is the spin label. So describes free electrons with 
chemical potential yit. 



Mx) , (3.2b) 



with m the fermion mass. We will find it useful to go to a Fourier representation with wavevectors k and fermionic 
Matsubara frequencies w„ = 2TrT{n+ 1/2), Q and a, {d + l)-vector notation, k — (k, a;„). Then Sq reads. 



k.a 



i^Ak) . (3.2b') 

2 m 



Here we have redefined S'o to include the first term on the right-hand side of Eq. ( 3.2a ). Sdis describes a static random 
potential, u(x), that couples to the fermionic number density, 

Sdis ^ - dx u{x)'^^pcr{x) tpa{x) , (3.2c) 
and Sint describes a spin-independent two-particle interaction, 

Sint ^ I dxidX2 V{xi - :ii2) ^paAxi)-4;a2{x2)i^a2{x2)-4^ai{xi) ■ (3. 2d) 

o-l,o"2 

The interaction potential w(x) will be discusse d fur ther below. 

For simplicity, the random potential in Eq. (3.2c) is taken to be delta-correlated with a Gaussian distribution. Its 
second moment is, 



1 



Mx)u(y)},^^ = J(x - y) , (3.3) 



where {. . .}^^^ denotes the disorder average. Here Td is the elastic-scattering mean-free time, and Np denotes the 
bare single-particle density of states (for both spin projecti ons ). More general random potential models can be used, 
but any differences between such generalizations and Eq. ( |3.3[ ) are irrelevant for the long-wavelength effects we are 
concerned with in this review. 

For disordered systems, single-particle momentum excitations are not long lived, but decay exponentially on a 
time scale given by the elastic mean-free time, Tci. The important physics on the longest length and time scales is 



We will denote fermionic frequencies by u!„, and bosonic ones by i}„ throughout. 
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controlled by two-particle excitations that are soft either because they are related to conserved quantities or because 
of a mechanism related to Goldstone's theorem. For a detailed discussion of the soft modes in a disordered, interacting 
electronic system we refer the reader to Ref. . The main conclusion is that of all the modes that appear in Sint , the 
dominant soft modes are those that involve fluctuations of either the particle number density n„, or the spin density 
n^, or density fluctuations, ric, in the particle-particle or Cooper channel. In Fourier space, these densities are given 
in terms of fermion fields by 

n„(g) = V77^IZ^^(fc)V'.(fc + 'Z) , (3.4a) 



n,(<z) = yrVF ^ V5.(fc)a..'V'.'(^ + 9) , (3.4b) 



fc 



'.cr.a 



with a ~ ((Tj., (Ty, cr^) the Pauli matrices, and 

nc(g) = v^E^T(fc)V'i(-fc + g) • (3.4c) 



This means that in Eq. (3. 2d) we should project onto these fluctuations. □ The net result is that Si„t can be split 



into three terms, each with a different interaction amplitude, 

Sint = slJf + S^J^ + S'|„| , (3.5) 
with the singlet, triplet, and Cooper channel terms simply describing interactions between the densities given in Eqs. 



(3.4), 



g 



C = ^E"^(9)-n.(-g) , (3.6b) 



(*) _ 

2 

9 



9 

HerePW, i = s,t, c, are singlet, triplet and Cooper channel interaction amplitudes, which are in turn given by angular 
averages of the interaction potential at specified momenta. 

Since the system contains quenched disorder, it is necessary to average the free energy or In Z . This is accomplished 
by means of the replica trick , which is based on the formal identity 

InZ = lim (Z" -l)/n . (3.7) 

Introducing n identical replicas of the system (with n an integer), labeled by the index a, and carrying out the disorder 
average, we obtain 

Z = {Z"}d.« = y7^[Vi",^1e^ . (3.8) 
We again separate S into free, disordered and interacting parts. 



®For short-ranged model interactions^ this can be done in a straightforward way [[in[ . A Coulomb interaction requires a 



somewhat more elaborate procedure ^Jj. However, the final result is given by Eqs. (3fB) in either case 
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S = J2(^o+Sls + Srnt) • (3.9a) 
Sq and J are given by Eqs. ( 3.2b1 ) and (^) and ( |3.6|) , respectively with ip ^ ^/j" and S*^-^ is given by, 



-^Ss = EEEE'^k,+k3,k.+k. c.(ki)c.(k2)viL.(k3)V'L,(k4) . (3.%) 

^ °' /3=1 {k} n,m cr,a' 

At the end of aU calculations, one takes the limit n —^ (assuming that both the analytic continuation from integer 
to real n and the limit exist) to obtain the desired averages for quenched disorder. 

B. Order parameter field theories 



Starting with the pioneering work of Landau |19|, it emerged over the years that an expedient way to describe any 
phase transition is to formulate an effective field theory, the LGW functional, entirely in term of an appropriate order 
parameter, and to bury the information about all other degrees of freedom in the vertices of that order parameter 
field theory. As we will see, this approach is more problematic for quantum phase transitions than for thermal 
ones. The reason is that the LGW approach hinges in part on an implicit assumption, namely that the order 
parameter fluctuations are the only soft modes in the system. Although this is not true for most of the quantum 
phase transitions we will be interested in, it turns out that the LGW approach nevertheless provides an economical 
way of solving the problem, provided one adequately deals with the ensuing complications. In the current subsection 
we therefore show how to derive LGW functionals of ferromagnetic, antiferromagnetic, and superconducting order 
parameters, respectively; for the time being without worrying about any additional soft modes. The question of an 
order parameter description of metal-insulator transitions is more intricate, and does require a thorough examination 
of the symmetry and soft mo de pr operties of the underlying fermionic field theory. Sec. Ill A, as a prerequisite. This 



analysis will be given in Sec. |III Q below. 

1. Magnetic order parameters 



It is well known that the model defined in Sec. [II A contains a ferromagnetic phase. Large values of the coupling 
constant P'-*-' and small values of T favor ferromagnetism, while small values of F*^*) and large values of T favor 
paramagnetism. The appropriate order parameter for the transition between these two phases is the magnetization 



m or the expectation value of the spin density defined in Eq. ( 3.4b| ): m = (ns(x)). Technically, the derivation of a 



LGW functional is easily achieved by means of a trick due to Stratonovich and Hubbard [gO|, which is based on the 
identity 



/ 



D[M] e-(r'*'/2)/''- M(.).M(.)+r(')/d. M(.).n.(.) ^ ^^^^^ ^ (g^^g) 



Here M(a;) is a classical (i.e., complex number valued) auxiliary field that couples linearly to the spin density. 
The strategy is now to apply this Gaussian decoupling procedure, which is referred to as a Hubbard-Stratonovich 
decoupling or as 'uncompleting the square', to S^^f. while leaving all other terms in the action unchanged. We thus 
write the action as 

S = Sref + In / D[M] e-(r'"/2) / M(.).M(.)+r<*) / d. M(.).n.(.) ^ (3^^^^ 



where the Sref contains all pieces of the action other than S^J^. We will refer to the system described by Sref as the 
'reference ensemble'. 

By applying the Hubbard-Stratonovich tra nsfor mation explained above, and formally performing the fermionic 



integral, the replicated partition function, Eq. (3.S), can be written as 



Z = e-''^°'^ J Yl D[M"] exp(-$[M]) , (3.12a) 
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with Fq the noncritical part of the free energy. The LGW functional <l> reads Q 



$[M] = -^Yl / M"(x) • M"{x 



In 



exp 



dx M"(a;) 



(3.12b) 



ref 



Here (. . .)re/ denotes an average taken with respect to the reference ensemble. A useful interpretation of the second 
contribution to the LGW functional is obtained by realizing that it represents the free energy of the reference ensemble 
in the presence of an external magnetic field that is proportional to the order parameter field M. Also notice that 
the average of M is proportional to the magnetization m, by virtue of the linear coupling between M and n^. 

For clean systems, and a more restricted reference ensemble, the above derivation of an LGW functional for quantum 
ferromagnets is due to Hertz [Q. Our derivation [|l| has generalized his approach to include disorder, and electron 
interaction in multiple channels. The standard way to proceed with a discussion of the paramagnet-to-ferromagnet 
transition is a Landau expansion or expansion of the LGW functional $ in powers of the order parameter field M. 
The structure of this expansion is 

oo 

1=2 

= \Y.I Xf;(xi,X2)M,"(xi)M,"(x2) 
a 

+ dxldX2dX3 X^f^{xi,X2,X3)M:{xi)M^{x2)M^{x:i) 

a 

- ^ E y dx2 dx3 dxi xillf{xi,X2,X3,X4) M^ixi) M^{x2)M^{x^) Aff (X4) + 0{M^) 



(3.12c) 



Q,/3 ■ 



The vertex functions X*^'^ are given in terms of the spin density correlations of the reference ensemble, and can be 
calculated within any theory of interacting disordered electrons. This we will come back to in Sec. [V below, where 
we will see that this Landau expansion is ill-behaved in the sense that the vertex functions cannot be localized in 
space and time. 

For an an tiferro magnet the interesting fluctuations are those of the staggered magnetization. Our simple continuum 
model. Sec. [HA, actually has no antiferromagnetic phase. However, order parameter theories for antiferromagnets 
can be deriv ed fro m more complicated microscopic theories The net result is again a LGW functional of the 

form of Eq. ( [3.12c ), with the fluctuating magnetization M replaced by a different vector order parameter N, whose 
average is proportional to the staggered magnetization. The structure of the ve rtex functions in this case is different 
from the ferromagnetic one. We will model these vertex functions in Sec. IV B. 



2. Superconducting order parameter 



For superconductors, the procedure is completely analogous to that in the ferromagnetic case, except that one uses 

(c) 

a Hubbard-Stratonovich decoupling on the Cooper channel interaction, Sl^^^. We thus write 

S — Sref + S-„j , (3.13) 



Fq is noncritical as long as the reference ensemble does not undergo a phase transition of its own. This is not an entirely 
trivial point, as perturbation theory within the reference ensemble will in general create a nonvanishing spin-triplet interaction, 
even though there was none in the bare system. For all of our order parameter theories we assum e that the reference ensemble 
is far from any phase transitions, since else the mode separation that is implicit in Eq. ( |j.ll[ ) breaks down, and nothing is 
gained by the Hubbard-Stratonovich transformation. 

^'^We denote the microscopic fermionic action by S, LGW functionals that depend only on the classical order parameter field 
by $, and effective fermion actions by A, respectively. 
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with 3^!^^ the rephcated generahzation of Eq. ( p. 6c ), and proceed as in Sec. IIIB 1, with the only difference being that 



we now have a complex valued order paramter that we denote by The LGW functional reads 



dq |^'"((j)p - In 



(3.14) 



ref 



The Landau expansion in powers of the order parameter is again analogous to the ferromagnetic case, Eq. (3.12c), 
except that due to gauge invariance only even powers of ^! appear, and due to the complex valuedness of ^' the vertex 
fmictions are scalars rather than tensors. This order parameter field theory will be discussed in Sec. ^ below. 



C. The nonlinear sigma— model 



Before we continue the discussion of the effective field theories derived in the previous subsection, it is useful 
to discuss the properties of the underlying fermionic field theory in some detail, in particular with respect to its 
symmetry properties and soft mode structure. This is very helpful for an understanding of the properties of the order 
parameter theories, and it is necessary in order to derive an order parameter theory for metal-insulator transitions. 
The symmetry analysis of the fermionic theory is very technical and involved, and the details can be found in Ref. 
[ psf . Here our goal is to explain the general structure of the theory, and the logic of the analysis, by exploiting a 
remarkable analogy between classical spin models for Heisenberg ferromagnets and disordered electrons that was first 
noted and used by Wegner p^ ]. 



1. Digression: The nonlinear sigma-model for classical Heisenberg ferromagnets 



Let us consider, as a model for a classical Heisenberg ferromagnet. 111 an 0{N) symmetric ^''-theory with a 
magnetic field h in the 1-direction. The action 



S[(j),h] = J 
determines the partition function 



r(0(x))'+c(V0(x))' 



u dx. 



(<?(x) • <?(x) 



dx 01 (x) 



Z[h] 



D[ 



'S[it>,h] 



(3.15a) 



(3.15b) 



and r, c, and u are real-valued coefficients that span the parameter space of the theory. The following discussion of 
this action, along standard lines p^], will turn out to be very useful as an analogy for the fermionic theory we are 
interested in. 



The crucial property of the action, Eq. (3.15a), that we want to exploit is its symmetry: For zero external magnetic 
field, /i = 0, the action is invariant under 0{N) rotations of the field 0. Let us consider a particular rotation, between 
01 and 0j; with 2 < i < N , through an infinitesimal angle 0, 



Pi ■ 



We also introduce a source J(x) for the 0-field, and consider the generating functional 

Z[h,j']^ I D\^]e~s[$.h]+Jd^J{^)-$i^) 



(3.16) 



(3.17) 



If we change integration variables from to 0', then the Jacobian of that transformation is equal to one. To linear 
order in the infinitesimal group parameter 9 we thus obtain 



'^■^The fact that classical ferromagnets are useful here for pedagogical purposes, while in Sees. [II Bl and IV A we discuss 
quantum ferromagnets, is accidental. 
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This source can be thought of as a spatial modulation of the magnetic field h. 
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/• 



D\$] ( -(5S'[(^,/i] + / dx J(x) • (5(^(x) ) e-^I-?'''! , (3.18) 



where 5S = S[(j)'] — S[(j)]. By differentiating with respect to Ji(x) and setting J = after the differentiation we obtain 
N — 1 Ward identities, 

(J5-0.(x)) + ((50,(x)) =0 , (i = 2,...,iV) , (3.19) 
where the brackets denote an average with respect to the action S. The only part of 5' that is not invariant under 



the rotation is the magnetic field term. The Ward identities, Eq. ( 3.19 ), thus take the form of a relation between 
the transverse (with respect to the magnetic field) two-point correlation functions and the longitudinal one-point 
function, 

hjdy (0,(x)0,(y)) = (0i(x)> . (3.20a) 

Since the left-hand side of this identity is h times the homogeneous transverse susceptibility xt , while the right-hand 
side is equal to the magnetization m, this can be rewritten as 

Xt = m/h . (3.20b) 

We thus obtain the well-known result that the transverse zero-field susceptibility diverges everywhere in the ordered 
phase, i.e. whenever m > 0. This is a particular example of Goldstone's theorem: Since there is an (N — l)-parameter 
continuous symmetry {0{N) in our case) that is spontaneously broken (m = {4>i)h^o 7^ in the ordered phase), there 
are — 1 massless excitations (the transverse (p-cj) correlation functions in the limit h ^ 0). 

In ordinary perturbation theory the masslessness of the — 1 Goldstone modes is not manifest, but rather has to 
be established by means of explicit calculations order by order j2^]. This is clearly undesirable, and raises the question 
of whether one can reformulate the theory so that this qualitative feature is explicitly displayed. This can indeed be 
done, and the result takes the form of a so-called nonlinear sigma-model |2^]. The basic idea is the realization that 
the massless Goldstone modes correspmid to purely transverse fluctuations of the vector fields, as the above derivation 
of the Ward identity demonstrates. Ej Let us therefore decompose the vector field into its modulus p and a unit 
vector field 0, 

0(x) = p(x) 0(x) , 02(x) = l , (3.21) 

Let V' be a particular fixed unit A'^-vector. Then any (/) can be generated from by means of an 0(Af)-rotation. The 
set of the is thus isomorphic to 0{N) modulo the subgroup that leaves ip invariant, which is 0{N — 1). The (j) 
therefore provide a representation of the homogeneous space 0{N)/0{N — 1). In terms of p and (j) the action reads, 

S[p,(f>]= J dx cp^{x) (\/${x)y + c{\/p{x)f + rp^{x) +u J dxp^{x)~h J dxp{x)^i{x) , (3.22a) 

and switching from the functional integration variables (/> to (p, 0) leads to a Jacobian or invariant measure 

I[p]=l[p''-\^) ■ (3.22b) 

X 

The important point is that in this formulation of the theory, the field (j) appears only in conjunction with two gradient 
operators. represents the — 1 soft Goldstone modes of the problem, while p represents the massive modes. Now 
we parametrize 0, 

0(x)-(a(x),7r(x)) , (3.23a) 

where 



Notice that the 0{N) transformation leaves the modulus of the vector field invariant. 
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a(x) = v/l-7?2(x) 



(3.23b) 



We split off the expectation value of the massive p-field, p(x) = i? + Ap(x), with R = (p(x)), and expand in powers 
of 7? and Ap. Rescahng the couphng constants with appropriate powers of R, the action can be written 



5[p, tt] = S_ 



NLaM 



[tt] + AS'[p,7r] 



Here 



SNLaAiln] = j I (V7r(x))' + (Va(x))' - h / dx(7(x) 



(3.24a) 



(3.24b) 



is the action of the 0{N)/0{N — 1) nonhnear sigma-model, which derives its name from the notation used in Eqs. 
(E23|), and 



AS'[p,7r] j {Ap{x)f+c J dx (VAp(x)) V 4i?u J dx (Ap(x))^ + u J dx (Ap(x))^ 



+ O Apa,Ap(\/(f> 



(3.24c) 



contains the corrections to it, which are all massive. If we neglect all fluctuations of the massive p-field, then we are 
left with the 0{N)/0{N — 1) nonlinear sigma-model in the usual parametrization. 

The separation of soft and massive modes is now complete, and will be preserved order by order in perturbation 
theory. We w ill exp lain in the next subsection how to perform an analogous separation of modes for the fermionic 
action of Sec. Ill A , Before we do that, however, let us briefly discuss a way to rewrite Eq. ( ^.24b| ) that will be very 
useful later. We enforce the constraint cjp' = 1 explicitly by means of a Langrange multiplier field A, 



Snl 



, A] = - / dx 



V0(x) +A(x) U2(x)-1 



(3.25a) 



Usually one integrates out A, which eliminates a in terms of tt, and expands in powers of tt . This leads to a (2 + e)- 
expansion for the Heisenberg transition. Alternatively, however, one can integrate out the tt fields to obtain an 
effective action Aeff in terms of a and A, 



Aeff [cr. A] 



dn. 



(V<7(x))2 + (a2(x) - 1) A(x)] + i (A^ - 1) Tr In (-V^ + A(x)) 



(3.25b) 



Now look for a homogeneous saddle-point solution. Denoting the saddle-point values of the fields again by a and A, 
respectively, we find 



Act 



^2^i_(iv-i)tly — ^ 

^ ^ V" ^ p2 + A 



(3.26) 



These saddle-point equations correctly reproduce the exact critical behavior for d > 4, which is mean-field like. They 
also yield the exact critical behavior for all d > 2 in the limit N ^ oo (with tN held fixed) (pij . 



2. Symmetry properties of the fermion model 



It is possible to rewrite the fermionic action of Sec. Ill A in complete analogy to the treatment of the much simpler 
(/)^-theory in the previous subsection. The first question that arises in this context is which symmetry group takes the 
place of 0{N). For this purpose, it is convenient to rewrite the fermionic theory in terms of four-component bispinors 
rj and 77+ that are defined by p7|] 



Ci(x) 
Ci(x) 



(3.27a) 
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and 



where 



with c the charge-conjugation matrix 



S2 
32 



(3.27b) 



(3.27c) 



(3.27d) 



The four degrees of freedom represented by the bispinor are the particle-hole or number density degrees of freedom, 
and the two spin degrees of freedom. We have also defined a basis in spin-quaternion space as Tr lE" Si (r, i = 0, 1, 2, 3), 



with To — So the 2x2 identity matrix, and Tj 



-iuj (j = 1,2,3), with fjj the Pauli matrices. In this basis. 



the channels r = 0, 3 and r — 1,2 describe the particle-hole and particle-particle degrees of freedom, and the channels 
1 = and i = 1, 2, 3 describe the spin-singlet and spin-triplet, respectively. In the space of bispinors we further define 
a scalar product, 



(3.28) 



Now let rj be transformed by means of an operator T: rj Tr]. By rewriting the action in terms of the bispinors, it 
is easy to see that the free fermion action Sq, Eq. ( |3.2b' ), is invariant under any transformations of the bispinors that 
leave the metric [rj,?]) invariant, except for the part that is proportional to the external frequency, iuj. The disorder 
part of the action Sdis, Eq. ( |3.9b ), is also invariant under such transformations. Remembering that 77"*" is related to rj 
by means of the charge conjugation matrix C, Eq. (3.27c), and using C"^ = , we find that in order to leave (77, rj) 
invariant, T must obey 



C 



(3.29) 



For a system with 2N frequency labels {N positive ones, including 0, and N negative ones), and n replicas, Eq. ( 3.29| ) 
defines a representation of the symplectic group Sp(8A^n,C) over the complex numbers C |^. Sp(8iVn,C) thus plays 
a role analogous to that of 0{N) in t^^-theoryj-and the external frequency in 5*0 can be thought of as being analogous 
to the external magnetic field in the latter. Ej The electron-electron interaction, which also breaks the symmetry, 
turns out to be very similar to iu in its effect, see below. 

The next question is whether there is a phase where the symplectic symmetry at zero frequency is spontaneously 
broken, and if so, what plays the role of the order parameter and the Goldstone modes. To see this it is convenient 
to define a Grassmannian matrix field -B12, with 1 = (rii,cki), etc, by 



Bi2(x) =r7+(x)®r72(x) , 
or, in Fourier space and with all components written out, 

.■'Bl(q)=E^'?nW^-77^,(k + q) 



(3.30a) 



(3.30b) 



where ^ry denotes the elements of t], and ^77 those of 77+. It is often convenient to expand the B into the spin-quaternion 
basis defined above. 



3 3 



B12 = E E 'rBl2 {Tr ® Si) 



(3.30c) 



r=0 i=0 



This is true even though iu couples to a term that is bilinear in or 77, while h in Sec. [II C 1 couples linearly to 
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The set of B is isomorphic to a set of classical, i.e. complex number valued, matrix fields. One can thus introduce 
a classical matrix field Q, and constrain _B to Q by means of a Langrange multiplier field, A. The fermionic degrees 
of freedom can then be integrated out, and one obtains a theory in terms of the matrix fields Q and A. In terms of 
these objects, the density of states, N, as a function of energy or frequency uj measured from the Fermi surface is 



iV(ej.+c.) = -Re(OQ„„(x 



Similarly, the number density susceptibility Xm and the spin density susceptibility, Xs, are given by 



mi ,m2 r- 



/ mi - 



.,(q);(AQ) 



m2,rn2 + 



„(-q) 



^0,3 



(3.31) 



(3.32) 



with AQ = g - (Q), and x^°^ = Xn and x 



/(1,2,3) 



Xs- = 27rTn is a bosonic Matsubara frequency. 



Under transformations of the bispinors, the Q and A transform accordingly. One can then derive Ward identities in 
analogy to Eqs. ( p.l^ ) - (3.20a), with the vector field (j) replaced by the matrix field Q, and the source J also being a 
matrix. The detailed calculation shows that the role of the transverse fields (* = 2, . . . , A'') is played by the matrix 
elements Qnm with nm < 0, i.e. by products of fermion fields whose frequency indices have different signs, while the 
analogs of the longitudinal field are matrix elements with nm > 0. The Goldstone mode equation analogous to Eq. 
( ^.20b|) takes the form. 



(k)^Q::;,(-k) 



k=0 



TrN{eF) 
16|fi„i_„2 



(nin2 < 0) 



(3.33) 



Here N^ep) is the exact density of states at the Fermi level. The salient point is that, as long as N{ep) > 0, the Q-Q 
correlation function at zero m omen tum diverges like l/|ri„j_„2|. We have therefore identified oQnin2 ^1*^2 < as 
a soft mode. As written, Eq. (3.33) is valid only for noninteracting electrons. However, an analysis of the inter action 
term shows that it does not spoil the property of the Q~Q correlation function on the left-hand side of Eq. ( 3.33| ) 
being massless as long as the density of s tates at the Fermi level is nonzero. Further, one can invoke additional 
symmetries of the action to show that Eq. (3.33) also holds for IQ with i, r 7^ 0. The soft modes in the particle-hole 
channel (r = 0,3) and particle-particle channel (r = 1,2) are often called 'diffusons', and 'cooperons', respectively. 
They are all soft in the absence of physical processes that break the additional sy mmet ries which link the general * Q 
to the qQ. Breaking these symmetries reduces the number of soft modes, see Sec. VIE below. 

It is important to realize that the above Q-Q correlation function corresponds to a very general four-fermion 
correlation, with no restrictions on the frequencies ni and 712 other than that they must have opposite signs. It 
is thus not related to the fermion number density, or some other conserved quantity, although the dens ity-d ensity 
correlation function can be obtained as a linear combination of these more general propagators, see Eq. ( 3.32| ). The 
physical reason for the softness of these modes is thus not a conservation law, but rather a spontaneously broken 
symmetry, viz. the rotation symmetry in -0 or 77-space between fermion fields with positive and negative frequency 
indices, respectively, or the symmetry between retarded and advanced degrees of freedom. This symmetry is broken 
whenever one has a nonzero density of states at the Fermi level (which is just the difference between the spectra of 
the retarded and advanced Green functions, respectively), and the soft Q excitations are the corresponding Goldstone 
modes. They are analogous to the transverse spin fluctuations in the classical Heisenberg model. 



3. Separation of soft and massive modes, and the nonlinear sigma-model for fermions 

We now know that the correlation functions of the Qnm with nm < are soft, while those with nm > are massive. 
Our next goal is to separate these degrees of freedom in such a way that the soft modes remain manifestly so ft to al l 



orders in perturbation theory, in analogy to the treatment of the 0(A'^)-symmetric Heisenberg model in Sec. IIIC 1 
Such a separation was first achieved by Schafer and Wegner for non-interacting electrons. The generalization of 
their treatment for the case of interacting electron was given in Rcf. ||l^. Here we sketch the logic of the reasoning, 
for details we refer to that reference. 

The matrices Q under consideration are complex 8Nn x 8Nn matrices, or, alternatively, quaternion-valued ANn x 
4Nn matrices. However, all of their matrix elements are not independent; the definition of the Q in terms of the 
fermionic bispinors implies symmetry properties that reduce the number of independent matrix elements. Using these 
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symmetry properties, |Ojrie can show that the Q can be diagonahzed by means of elements of the unitary symplectic 
group USp(8iVn,C). M The most general Q can thus be written 



Q = S D S^^ , (3.34) 

where D is diagonal, and S € USp(8A^n, C). 

However, diagonalization is more than we want. Since we know that the Qnm with nm < are soft, while those 
with nm > are massive, we are interested in generating the most general Q from a matrix P that is block-diagonal 
in Matsubara frequency space, 

^ = (T P<) ' (^-^^^ 

where and P< are matrices with elements Pnm where n,m > and n,m < 0, respectively. This can easily be 
achieved. Due to the symmetry properties of and P^, the most general P can be obtained from D by an element 
U of USp(4A^n,C) x \JSp{4:Nn,C). The most general Q can therefore be written 

Q = S P S-^ , (3.36) 

with S = SU^^. The set of transformations S is the set of all coscts of \JSp{8Nn,C) with respect to lJSp{4Nn,C) x 
USp(4iVn,C), i.e. the S form the homogeneous space USp(8Arn,C)/USp(4iVn,C) x USp(47Vn,C). 

This achieves the desired separation of our degrees of freedom into soft and massive ones. The massive de- 
grees of freedom are represented by the matrix P, while the soft ones are represented by the transformations 
S £ USp(8iVn,C)/USp(4Afn,C) x USp(47Vn, C). The analogy with the 0(iV)-Heisenberg model is now obvious: 
The unitary-symplectic coset space, identified above as the space whose elements comprise the soft modes of the 
theory, is a matrix generalization of the homogeneous space 0{N)/0{N — 1) x 0(1), which represents the soft modes 
(j) in the 0{N) vector model. 

In order to formulate the field theory in terms of the soft and massive modes, one also needs the invariant measure 
I[P], or the Jacobian of the transformation from the Q to the P and the S, defined by 

J D[Q] ... = J D[P] I[P] J D[S] ... . (3.37) 

I[P] has been constructed explicitly p^ . However, it is not necessary to know its explicit form to derive the nonlinear 
sigma-model. 



The further development of the theory is analogous to the derivation of Eqs. (3.24), although the technical details, 
are more intricate [|l^. One splits the expectation value (P) off the massive field P, and defines a new matrix field tZi 

Q(x)=5(x)(P)5-i(x) , (3.38) 

Q has the properties = const, and trQ = 0. The action can then be written as a generalized matrix nonlinear 
sigma-model that contains only the soft modes Q, and corrections that are all massive. The nonlinear sigma-model 
part of the action reads 

ANLaM = ^ y dx tr (vQ(x))^ + 2H J dx tr (f7Q(x)) + Ant[Q] , (3.39a) 

with the matrix field Q — Q — _(P), and Amt the interaction part of the action, Eq. ( |3.5D , replicated and rewritten 
in terms of the Q-matrices. c3 The explicit derivation yields the coupling constants as C? = B/ttcto, with (Jq the 
conductivity in self-consistent Born approximation, and H = irNp/S which can be interpreted as the bare quasi- 
particle or specific heat density of states j3^,^,^. G serves as a measure of the disorder in the system. 51 is a 
diagonal matrix whose diagonal elements are the fermionic Matsubara frequencies a;„. Q is subject to the constraints. 



^^This group is defined as the intersection of a unitary and a symplectic group, USp(2n,C) = U(2n,C) nSp(2n,C) [ggl- It can 
be represented by unitary matrices that are also symplectic. 
'^^Strictly speaking, Q contains only the real part of (P), and we also drop a normalization factor. 

'^^One can write either Q and Q in the action, Eq. (3. 39a). The only resulting difference is an uninteresting constant 
contribution to AnlitM- 
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Q2(x) = 1 , Qt = Q , trQ(x)=0 



(3.39b) 



The effective action given by Eqs. ( 3.39 ) is the generahzed nonhnear sigma-model that was first proposed by 
Finkel'stein | |33| , p^ as a model for disordered interacting electrons, and whose properties have been studied in consid- 
erable detail, in particular with respect to the metal-insulator transition that is described by it l35| , |3^ , ^ , ^ , ^ . We 
will review the properties of this theory in Sec. below. 



4- Order parameter field theory for metal-insulator transitions 



From Eq. ( 3.33| ) we know that the order parameter that is related to the existence of Goldstone modes in the 
Fermi system is the density of states at the Fermi level, N{eF)- It is known that at the metal-insulator transition 
of interacting electrons, N^ep) is critical [Q. It is nonzero in the metallic phase, and vanishes as the transition is 
approached. This suggests that N{€p) can be considered an order parameter for this type of transition, and that one 
should try to derive an order parameter field theory for t he me tal-insulator transition in analogy to those for the 
magnetic and superconducting transitions discussed in Sec. IIIB . One way of doing this is as follows |37y3§f| . 

In terms of the Q-matri x, N( ep) is given as the expectation value of the diagonal matrix element Qnn in the limit 
of zero frequency, see Eq. ( ^.31 ). A pure order parameter description would therefore integrate out all modes except 
for Qn=o,n=o- Since the frequencies are coupled, however, it is adavantageous to keep all of the modes that are massive 
in the ordered (i.e. metallic) phase, and to integrate out the soft modes, i.e. the Qnm with nm < 0. This can be done 



(3.40) 



in analogy to Eqs. ( p.25| ). We separate Q into blocks in frequency space, 

, if n > 0, TO > or ri < 0, m < 
, if n > 0, TO < , 

if 71 < 0, TO > . 




and enforce the constraint = 1, Eq. ( 3.39b ), by means of a matrix Lagrange multiplier field A, 
n S[Q^^) -l]= j D[K] exp{-^ j dx tr (A(x)[Q2(x) - 1])} , 



(3.41) 



where the factor of 1/26* has been inserted for convenience. By decomposing A int o blocks like Q one sees that the 
elements A„m with nm > 0, together with the tracelessness condition, Eq. (3.39b), are sufficient for enforcing the 
constraint = I. We can therefore restrict ourselves to A„m with nm > 0. The action, Eq. ( p.39a ), can then be 
written 



A(x)[Q2(x) - 1] + (vQ(x)y 



2H / rfx tr 



(r!Q(x)) +A^.. 



(3.42) 



The soft modes q can now be integrated out to obtain an effective action entirely in terms of Q and A. We obtain 



Aeff^ln J D[q,q'<] CXp(ANLaM[Q,M 



A 



NLcrM 



[Q, A] - - Tr lnM[A] 



(3.43) 



Here M is a complicated matrix that depends linearly on A. By constructing a saddle-p oint solution of Aeff one can 
now construct a mean-field theory of the metal-insulator transition in analogy to Eqs. ( 3.26| ). This will be discussed 
in Sec. |VIB2| below. 



IV. MAGNETIC TRANSITIONS AT ZERO TEMPERATURE 



A. Itinerant ferromagnets 

Perhaps the most obvious example of a quantum phase transition is the paramagnet-to-ferromagnet transition of 
itinerant electrons at T = as a function of the exchange interaction between the electron spins. Early theoretical 
work on this transition suggested that the critical behavior in the physical dimensions d — 2 and d — 3 was mean-field 
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like, as it is in the thermal ferromagnetic transition in d > 4. As discussed in Sec. |, the reason for this suggestion was 
that the phase transition in a quantum mechanical system in d dimensions is related to a classical phase transition 
in d + z dimensions. Formally it was shown Q that z — 3 and z = 4 in clean and disordered itinerant quantum 
ferromagnets, respectively. This appeared to reduce the upper critical dimension djt, above which fluctuations are 
unimportant and simple mean-field theory yields the correct critical behavior, from = 4 in the classical case to 
d^ — 1 and (i+ = in the clean and disordered quantum cases, respectively. If this were true, then this QPT would 
be uninteresting from a critical phenomena point of view. 

It is now known that this is not the case ||2l| , |39| , ^0| . In the early work there was an implicit assumption that the 
only modes or excitations that are important for describing_the phase transition were the critical modes. For the finite 
temperature phase transitions this is generally correct. However, at zero temperature the statics and dynamics 
are coupled together, which tends to increases the number of soft modes. Further, at T = there are soft excitations 
that exist in the entire metallic phase due to a spontaneously broken symmetry that is unrelated to magnetism. These 
particle-hole excitations ('diffusons') and particle-particle jescitations ('cooperons') were discussed in Sec. III. These 



modes cause what are known as weak localization effects, l£j and they are soft due to their diffusive nature. It turns 
out that these additional soft modes couple to the critical modes and influence the critical behavior. 



1. Disordered ferromagnets 



We now proce ed wi th the discussion of the LGW functional derived in Sec. IIIB , along the lines of Ref. A s 
discussed in Sec. IIIB , formally expanding the LGW functional ^[M], Eq. ( 3.12b|) , m powers of M yields Eq. ( |3.12c| ). 
Of particular interest to us is the quadratic or Gaussian term, which in Fourier space is given by 



(4.1) 



q,f2„ a 



Here we have scaled M with 1/T^*\ Xs is the Fourier transform of the dynamical spin susceptibility in the reference 
ensemble. 



Xsixi - X2) = {nlJxi)nl^{x2))^^^ , 



(4.2a) 



where rig, a (^ = 1, 2, 3) is one component of the spin density vector n^, Eq. ( 3.4b| ). Spin density conservation implies 
that at small frequency and wavenumber, the Fourier transform of Xs has a diffusive structure 



Xs{q,^n) = xo(q) 



Dq^ 



(4.2b) 



where D is the spin diffusion coefficient and Xo(q) is the static spin susceptibility, both in-the reference ensemble. In 
order to reach criticality, the frequency must be taken to zero before the wavenumber. Ell In the critical limit, we 
can thus expand. 



Xs(q, aO = Xo(q) [l-\n,,\/Dq^ + ...] 



(4.2c) 



It can be easily verified that the corrections to the leading terms in Eq. (4.2c) are irrelevant for the critical behavior. 



However, complicated soft mode structures can occur in classical systems as well, for instance in liquid crystals. See Ref. 
[ pl| for an example. 

The term 'weak localization' is loosely defined, and often misunderstood. We use it to refer to the nonanalytic behavior of 
electronic correlation functions in the limit of zero momentum and/or frequency that is induced by quenched disorder, or by a 
combination of interactions and quenched disorder, and that occurs even if the system is far from any kind of phase transition. 
The physical reason for these nonanalvticities is the difltusive motion of the electrons in the presence of quenched disorder. For 
a detailed discussion see, e.g., Refs. |^^, and p^ . 

^^This can be seen as follows. Since the magnetization is conserved, ordering on a length scale L requires some spin density to 
be transported over that length, which takes a time t ~ /D, with D the spin diffusion coefficient. Now suppose the coherence 
length is ^, and we look at the system at a momentum scale g or a length scale L ~ 1/g < ^. Because of the time it takes the 
system to order over that scale, the condition for criticality is < Min(_Dt, ^^). In particular, one must have < Dt, or 
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Next comes the crucial step of calculating Xo(q) in the reference ensemble. For a long time it was assumed that 
Xo(q) was an analytic function of q^. This might seem plausible, given that the reference ensemble describes a physical 
system that is far away from any critical point. However, it is now known that this reasoning is incorrect, and that 
any itinerant electron system at zero temperature has long-range correlations everywhere in its phase diagram, even 
far away from any critical point Q]. These long-range correlations are due to the extra soft modes that generically 
exist in quantum systems. In disordered systems, these effects are known collectively as weak localization effects. For 
small values of |q|, and for d > 2, the leading behavior of xo(q) is 



xo(q) 



Co - Cd-2|q 



d-2 



C2q 



(4.3) 



where the Ci are positive constants. Notice that the susceptibility decreases with increasing wavenumber. The 
negative sign of the jq]'^"^ term, whose prefactor Cd-2 vanishes for noninteracting electrons, is due to the fact that 
the interaction increases xo, and its effect gets weaker at larger values of |q|. For d < 2, the electrons are localized 
and a different theory is nee ded. 

Using Eqs. (4^) and (O), the Gaussian part of the LGW functional can be written more explicitly. 



$2(M) 



E2:[^o+q' 



d-2 



|f}„|/q21|M"(q,f], 



q,a„ 



where we have omitted the prefactors of the various terms, and 



1 



Xo(q = 0) 



(4.4a) 



(4.4b) 



is the bare distance from the critical point. We recog nize E g. (4.4b) as a generalized Stoner criterion for ferromag- 
netism. Physically, the most important term in Eq. ( 4.4a ) is the |q|'^^'^, which means that in real space the spin 
density fluctuations interact through a potential that decays with distance r as r^-^^,. It is well known from the 
theory of classical phase transitions that long-range interactions suppress fluctuation effects, and the critical behavior 
in systems with such interactions can therefore be determined exactly |4^ . For example, using renormalization group 
methods it can be established that all terms higher than I = 2 in the Landau expansion, Eq. ( [3.12c ), are renormal- 
ization group irrelevant for c? > 2, so that the upper critical dimension for fluctuations is d^ =2. From this it follows 
that many of the critical exponents can be determined exactly from $2. 
To see this we define the order parameter correlation function by 



G(q,f^„) = (|A4(q,a^)P) , 
with Ma an arbitrary component of the vector field M. For small |q| and f2„ we find 



G(q, r!„) = 



1 



d-2 



(4.5a) 



(4.5b) 



Here we have again omitted all prefactors of the various terms in the denominator in Eq. ( 4.5b| ), since they are 
irreleva nt fo r our purposes. In this equation we have also anticipated that irrelevant variables will renormalize the to 
in Eq. ( 4.4a| ) to t, the actual distance from the critical point, in Eq. ( 4.51 ). For q = Sl„ = 0, the correlation function 
G determines the magnetic susceptibility Xm ^ G(q = 0, 17 = 0) in zero field. Hence we have Xm{t) ~ = , 
where the last relation defines the critical exponent 7. In our case. 



7 



(4.6) 



which is valid for all d > 2. 7 thus has its usual mean-field value. However, for non-zero wavenumbers the anomalous 
|q|"'~^ term dominates the usual q^ dependence for all d < A. By inspection, the correlation length exponent, defined 
by ^ ^ t""^, is given by 



l/{d- 
1/2 



for 2 < d < 4 
for d > 4 . 



(4.7) 



Note that v > 2/d, as it must be according to the discussion in Sec. ||. The wavenumber dependence of G at t = is 
characterized by the exponent r/, which is defined as G(q, fin — 0) ^ |q| 



'2+rt 



From Eq. (4.5b) we obtain 



4-d , for 2 < d < 4 
, for d > 4 . 
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Finally, the dynamical scaling exponent z is read off the relation ^ ^ jt to be 



d , for 2 < d < 4 , , . 

4 , for d > 4 . ^ > 

Note that all of these exponents 'lock into' their mean-field values at a dimensionality fi++ — 4. 

To determine the critical exponents /? and 8 we need the equation of state in the ordered phase. An important 
point is that although the quartic {i = 4) term is irrelevant in the technical sense of the renormalization group, it 
is nevertheless needed to determine the critical behavior of the magnetization. The reason is that the coefficient 
of the M'*-term is an exa mple of a dangerous irrelevant variable. Further, the same mechanism that causes the 

lead s to a divergence of that coefficient in the long-wavelength, low-frequency limit, i.e. 



nonanalyticity in Eq. (4.3 



the function X^'^'> in Eq. ( 3.12c ) cannot be localized in space and time. The terms with Z > 4 in Eq. ( 3.12c ) are even 
more singular in this limit. Despite this unpleasant behavior of the field theory, it is easy to see by power counting 
that all of these terms are irrelevant with respect to the Gaussian fixed point. Taking these singularities into account 
in constructing the equation of state gives for d > 2, 

tm + m'^/'^ + rr? ^ h , (4.10) 

where we again have left off the irrelevant prefactors. Here m is the modulus of the magnetization vector, and h is the 
external magnetic field. Notice the term m''/^, which occurs in addition to what is otherwise an ordinary mean-field 
equation of state. It occurs because of the singularities mentioned above. For d < 6 it dominates the nr' term, and 
hence determines the critical exponents /3 and 6. The zero-field m agne tization near criticality behaves by definition 



of the exponent /3 as m{t, /i = 0) ~ . The equation of state, Eq. (4.10), gives 



Similarly, the exponent 5, defined by mit — 0,h) ^ h^/^ , is 

^^rd/2 , for2<d<6 , 
1^ 3 , for d > 6 . 

Note that these relations imply yet another upper critical dimension, rfj^^^ = 6, defined as the dimension where 13 
and 5 'lock into' their mean-field values of 1/2 and 3, respectively. 

The critical behavior of the specific heat, cy, has also been calculated. It is most convenient to discuss the specific 
heat coefficient, 7y — lim^^o cy/T, which in a Fermi liquid would be a constant. Its behavior at criticality, i = 0, is 
adequately represented by the following integral, 

Jo T + + + ^/"jq^ 

Remarkably, in zero magnetic field, 7y diverges logarithmically as T ^ for all dimensions 2 < d < 4. From a scaling 
viewpoint this is a consequence of the dynamical scaling exponent z bein g exactly equal to the spatial dimensionality, 
d, in this dimensionality range. In terms of a scaling function, Eq. ( 4.12| ) implies 



7y(i,T,/i) = 6(4- fi)ln6 + FT,(i6i/'', . (4.13) 

Here 0(x) denotes the step function, h is an arbitrary scale factor, and is a scaling function that cannot be 
determined by general arguments alone. 

The complexity of the critical behavior displayed above is astonishing, given that it is controlled, after all, by 
a simple Gaussian fixed point. For instance, the appearance of several different 'upper critical dimensionalities' is 
suprising at first sight. As the scaling analysis makes clear, this is a direct consequence of the nonanalytic terms in 
the LGW functional, which in turn results from the 'additional' soft modes. Furthermore, the critical behavior of 
different quantities can be affected by different dangerous irrelevant variables. For instance, the critical behavior of the 
magnetization is governed by the coefficient of the (irrelevant) quartic term in the LGW functional, while the specific 
heat coefficient is not affected by this dangerous irrelevant variable. As a result, the Wilson ratio, W — {m/h)/^Vi 
diverges at the criticality, while general scaling predicts that it should be a universal finite number. 

Another interesting example for the breakdown of general scaling is provided by the temperature dependence of 
the magnetization m, and of the magnetic susceptibility Xm- A detailed analysis pf^jpOl reveals that these quantities 
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do not depend on the critical temperature scale, and therefore their temperature scaling behavior is not determined 
by the dynamical critical exponent but rather by the scale dimensions of subleading temperature scales. One finds 
for m and Xm the following homogeneity laws, 

m(t, r, H) = h-^^'^mitb^/", Tfe-^/^ Hh^P''') , (4.14a) 

Xm{t, T, H) - K''''xm{th^'\Tb^'\ Hb'^'l ■ (4.14b) 
Here the crossover exponent (j) turns out to be given not by vz, but rather by @ 

2v^2/{d-2) , for 2<d<y/l+l , 

d/2 , for V5 + 1 < d < 4 , ^^^n 

4/(d- 2) , for 4 < d < 6 , V • ^ 

2v = l , for d > 6 . 

The correct physical interpretation of the temperature dependence of m and Xm is not dynamical scaling, but rather 
the crossover from the quantum critical region to a regime whose behavior is dominated by the classical Gaussian 
fixed point [|§. 

2. Clean ferromagnets 

In the previous subsection, the problem of disordered ferromagnets was considered. The only point in that dicussion 
where the disorder was important was the diffusive dispersion relation of the 'extra' soft modes. This raises the question 
whether similar effects might exist in clean itinerant ferromagnets. The first question that arises in the context is 
what, if anything, will replace the q''~^-term in the static spin susceptibility, Eq. ([4.3|), in clean syst ems. To answer 



this, let us consider the perturbation theory for Xo(q)- What leads to the nonanalyticity in Eq. (O) is the coupling 



of two diffusive modes, which mathematically takes the form of a mode-mode coupling integral of the type 

f 1 1 

dk / du) — — , (4.15) 

7 w + k2tj+(q + k)2 ' ^ ^ 

with q the external momentum, and we have set the diffusion coefficient equal to unity. Renormalization group 
techniques have been used to show that this is indeed the leading small-q behavior of xo- 

What changes in a clean system? The soft modes are still the density and spin density fluctuations, and in addition 
more general particle-hole excitations. All of these have a linear dispersion relation, i.e., lo ^ |q|. One might thus 



expect Xo(ct) in a- clean system to have a mode-mode coupling contribution analogous to that given by Eq. (4.15), 
but with ballistic modes instead of diffusive ones, 

dk / duj ^— ^ . (4.16) 

7 tj+ |k| a;+ |q + k| 

In generic dimensions, expanding in |q| leads to jslj], 

Xo(q)- const + dd_i|q|''-i+d2|q|2 , (4.17) 



at T = 0. Here dd-i and d2 are constant prefactors. For d < 3, the nonanalytic term in Eq. (4.17) represents the 
leading q-dependence of xo- In d = 3, one finds a q^lnl/jql term, and in d > 3 the analytic contribution is the 
leading one. 



If dd-i < in Eq. (4.17), then the same arguments used in the previous section apply here. For d < 3 one obtains 
the critical exponents 

/5 = ^ l/(d- 1) , 77 = 3-d , S = z^d , 7 = 1 , (4.18) 



^Refs. |2^,^ contained a mistake in the result for 0, which was corrected in Ref. 
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and Eq. (4.13) is valid with 4 — d replaced by 3 — d. For d > 3 all of the exponents have their mean-field values, and 
in d = 3 there are logarithmic corrections to mean-field like critical behavior. 

As in the disordered case discussed in the previous subsection, the temperature scaling of the magnet iz ation i s 



complicated and violates general scaling. The homogeneity laws for m and Xm are given by Eqs. ( 4.14a ), ( 4.14b ) 
with a crossover exponent (j) given by ||3^ 

!/==l/(d-l) , forl<d<2 , 

d/2(d-l) , for2<d<3 , , 

3/(d+l) , for3<d<5 , ^ ' 

v=l/2 , ford>5 . 

These results for the quantum critical behavior of clean itinerant ferromagnets disagree with those of Ref . Q| . It 
should be emphasized that these discrepancies cannot be related to differences in the underlying models. If the model 
of Ref. is properly renormalized, then it shows the same critical behavior as discussed above, since all features of 
the action that are not included in the bare model will be generated by the renormalization group. This was shown 
in Ref. Jsot . It highlights the important point that a simple power counting analysis yields the right answer only if 
no qualitatively new terms are generated under renormalization. 



Detailed calculations have confirmed the existence of the |q|''~^-term in Eq. (4.17) |49|. However, to lowest order 
in perturbation theory the coefficient dfi^i of that term is positive. For d < 3 this implies that xo increases with 
increasing q like Iq]''"^. For any physical system for which this were the true asymptotic behavior at small q, 
this would have remarkable consequences for the zero-temperature phase transition from the paramagnetic to the 
ferromagnetic state as a function of the exchange coupling. One possibility is that the ground state of the system will 
not be ferromagnetic, irrespective of the strength of the spin triplet interaction, since the functional form of xo leads 
to the instability of any homogeneously magnetized ground state. Instead, with increasing interaction strength, the 
system would undergo a transition from a paramagnetic Fermi liquid to some other type of magnetically ordered state, 
most likely a spin density wave. While there seems to be no observational evidence for this, let us point out that in 
d = 3 the effect is only logarithmic, and would hence manifest itself only as a phase transition at exponentially small 
temperatures, and exponentially large length scales, that might well be unobservable. For d < 2, on the other hand, 
there is no long-range Heisenberg ferromagnetic order at finite temperatures, and the suggestion seems less exotic. 
Furthermore, any finite concentration of quenched impurities will reverse the sign of the leading nonanalyticity, and 
thus make a ferromagnetic ground state possible again. 

Another possibility is that the zero-temperature paramagnet-to-ferromagnet transition is of first order. It has 
been shown in Ref. p9| | that the nonanalyticity in xo {^) leads to a similar nonanalyticity in the magnetic equation of 
state, which takes the form 

tm-Vdm'^ + um^ = h . (4.20) 

with m the magnetization, h the external magnetic field, and u > a positive coefficient. If the soft mode mechanism 
discussed above is the only mechanism that leads to nonanalyticities, then the sign of the remaining coefficient Vd in 



Eq. (4.20) should be the same as that of dd-i in Eq. (4.17), i.e. Vd > 0. This would imply a first order transition for 



2 < d < 3. In this case the length scale that the scenario of the previous paragraph would have been attributed to a 
spin density wave would instead be related to the critical radius for nucleation at the first order phase transition. 
Finally, it is possible that terms of higher order in the interaction amplitude could change the sign of dd-i in Eq. 



(4.17) for interaction strengths sufficient to cause ferromagnetism. In this case one obtains a continuous quantum 



phase transition with critical exponents as given by Eqs. (4.18) and (4.19). Further work is clearly needed on the 
clean itinerant quantum ferromagnetism problem. 

B. Disordered antiferromagnets 

Quantum antiferromagnetism (AFM) has experienced a surge of interest in recent years, mainly in efforts to 
explain the magnetic properties of high-Tc materials. Much of the theoretical work has focused on one-dimensional 
spin chains, both with and without quenched disorder. In d = 3, the QPT in clean itinerant systems is mean-field 
like, since z — 2 and hence d + z is greater than the classical upper critical dimension djf = 4 for all d > 2. In this 
case Hertz's LGW theory does not break down, since the order parameter (the staggered magnetization) is not a 
homogeneous quantity. Therefore the additional soft modes that make the ferromagnetic problem so interesting do 
not couple to the order parameter fiuctuations, and the vertices in Eq. ( 3.12c[ ) can be localized in space and time. 



There has been relatively little work, on the other hand, on the disordered AFM problem in d > 1. Even though in 
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this case as well the soft noncritical modes of the previous section do not influence the AFM transition, the inclusion 
of quenched disorder leads to a very interesting and difficult problem. 

Following earlier work by Das Gupta and Ma , Bhatt and Lee |5l[] have studied some aspects of the interplay be- 
tween AFM and strong quenched disorder. They considered a model for the insulating phase of a doped semiconductor 
consisting of an ensemble of randomly distributed, AFM coupled Heisenberg spins with a very broad distribution of 
coupling constants, J. Using a numerical renormalization procedure they found that with decreasing temperature, an 
increasing number of spin pairs freeze into inert singlets, but that some unpaired spins remain and give an essentially 
free spin contribution, Xm ~ ^/T, to the magnetic susceptibility as T — > 0. The net result was a 'random-singlet' 
phase, with a sub-Curie power law T dependence of the magnetic susceptibility as T — > 0. The quantum nature of 
the spins thus prevents the classically expected long-range order of either AFM or spin glass type. 

Bhatt and Fisher |Q applied similar ideas to the disordered metallic phase. They argued that rare fluctuations in 
the random potential always lead to randomly distributed local moments to which the methods of Bhatt and Lee can 
be applied, but now in the metallic phase. Their conclusion was that the local moments cannot be quenched by either 
the Kondo effect or the RKKY interaction induced by the conduction electrons. This again leads to a random-singlet 
phase with a magnetic susceptibility that diverges as T — > 0, albeit slower than any power. 

These results seem to imply the unlikely conclusion that antiferromagnetic long-range order can never exist in the 
ground state of disordered systems. In d = 1 this conclusion is known to be wrong |]5^ . In higher dimensions it clearly 
warranted further investigation. In Rcf. [ p4[ it was established that in some parts of the phase diagram, long-range 
order can exist. The basic idea is that quantum fluctuations weaken the metallic random-singlet phase since they 
enhance the interaction of the isolated local moment electrons with their environment. These interactions in turn 
restore long-range order by suppressing the random-singlet phase that would otherwise pre-empt an AFM transition. 

This result was established by studying a model of an itinerant AFA/I-a/ith a spatially random Neel temperature, or 
mass term. This model can be derived from a fermionic description, 123 just like the ferromagnetic one was derived 



in Sec. III. The only difference is that the antiferromagnetic order parameter, N(x, r), is multiplied by a phase 
function that represents perfect AFM order. Because of this phase factor the singular soft modes responsible for the 
nonlocalities in the ferromagnetic LGW theory do not couple to the order parameter. The resulting replicated local 
field theory for N has an action 

5[N"] = i /dxdy f ^ dTdr'V] N"(x,T)r(x-y,T-r')N"(y,r') 
Jo 

+u / dx / dr V(N"(x,t) •N"(x,t))^ 
J Jo 

-A [ [ drdr' ^ (N"(x,t) •N"(x,r)) (N'^(x,t') •N''(x,r')) . (4.21) 

The Fourier transform of the bare two-point vertex function r(x, r) is 

r(q,f]„)-t + q2 + |r!„| . (4.22) 

u is the coupling constant for the usual quartic term, and A is the strength of the random mass term. Physically, u 
is a measure of the strength of quantum fluctuation effects, while A is a measure of disorder fluctuations. Note that 
w, A > 0, so the presence of disorder (i.e., A 7^ 0) has a destabilizing effect on the field theory. To initiate the analysis 



of Eq. (4.21) we put T — 0, and first consider the case A = 0. We again define the scale dimension of a length L 



to be \V\ — —1 (see Eq. (2.5)), and that of time to be [r] = — z, and look for a Gaussian fixed point where z — 2 
and 77 = 0. Power counting in d dimensions shows that \u\ = 4 — (d + z), so that u is irrelevant, and the Gaussian 
fixed point is stable, for all d > 2. The other exponents can then be readily determined; for instance, v = 1/2. In 
contrast, the term proportional to A carries an extra time integral, so with respect to the Gaussian fixed point we 
have [A] = 4 — d. Hence the disorder is relevant for d < 4, and the Gaussian fixed point is no longer stable in the 
presence of disorder. This instability of the Gaussian fixed point can also be inferred from the Harris criterion 

In order to see if there are other fixed points that might be stable, a one-loop renormalization group calculation 



for the model given by Eq. (4.21) was performed in Ref. ]54| . The results are shown in a schematic flow diagram on 
the critical surface in Fig. One finds indeed a stable, critical fixed point that we call the random fixed point. In 



^^Of course it is necessary to consider fermions on a lattice, as a continuum model does not display antiferromagnetism. 
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FIG. 2. Schematic flow diagram on the critical surface. The unstable Gaussian (G) and the stable random (R) fixed points 
are shown. The dotted line denotes the boundary of the basin of attraction for the random fixed point. (From Ref. [p4|.) 



Figure both the Gaussian and the random fixed points are indicated. The linear renormalization group eigenvalues 
on the critical surface are complex, so that the corrections to scaling at the transition corresponding to the random 
fixed point are oscillatory in nature. More importantly, the random fixed point has only a limited range of attraction. 
If the initial values of A and u are inside the basin of attraction for that fixed point, then the system undergoes a 
continuous phase transition, and standard techniques can be used to calculate the critical exponents. If, on the other 
hand, the initial values of A and u place the system outside of that basin of attraction, then one finds runaway fiows. 
The correct interpretation of this result is not obvious. One possible interpretation is a fluctuation induced first order 
phase transition. However, for a state with simple AFM order this possibility is not realized |E5|. Another obvious 
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FIG. 3. Schematic phase diagram, showing the antiferromagnetic (AFM) and random-singlet (RS) phases in the J — A 
plane, with J the effective AFM coupling constant, and A from Eq. (4.21). For A = 0, one has an AFM for J > Jc, and a 
Fermi liquid for J < Jc. The dependence of the phase boundary on the fluctuation parameter it is not shown. (From Ref. [M.) 



possibility is that for some parameter values, long-range order is not possible. The latter possibility is consistent with 
the phase diagram shown in Fig. ||, in which the random-singlet phase (no long-range order) and the AFM phase 
compete with one another. The basic idea is that although disorder destabilizes the field theory, and hence works 
against ordering, quantum fiuctuations as measured by u can counteract this and lead to long-range order. For this 
to happen, however, the initial values at u and A must not be too small since u is a renormalization group irrelevant 
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variable for d = 4 — e. In other regions of parameter space where no critical fixed point exists, the disorder wins out 
and scales to large values. The resulting field theory is unstable, but has local instanton or local magnetic moment 
solutions. Presumably these local solutions represent the random singlet phase. Qualitatively, the most striking 
prediction is that for small disorder, there is no long-range order, but that for intermediate values of A, long-range 
order does exist. For further discussions we refer to the original publication [M. 



V. SUPERCONDUCTOR-METAL TRANSITION AT ZERO TEMPERATURE 



The transition from a normal metal to a superconductor at T = provides an interesting contrast to the magnetic 
transitions discussed in Sec. [ p6[ . Again there are soft modes in addition to the order parameter fluctuations, 
in this case particle-hole excitations in the Cooper channel or 'cooperons', that are integrated out in deriving a 
LGW functional, and that lead to long-range effective interactions between the order parameter fluctuations. For 
superconductors, the additional soft modes are stronger than in the magnetic case, and they completely dominate the 
physics. As a result, the quantum critical behavior of superconductors is BCS-like. There are, however, very strong 
corrections to scaling that lead to a broad transition region as d ^ 2. 

Let us concentrate on the realistic case of a superconductor with a nonzero density of elastic scatterers, so that the 
normal state conductivity is finite. The transition can then be triggered by varying either the attractive electron- 
electron interaction, or the disorder. We proceed by performing a Landau expansion or expansion of the LGW 



fimctional $, Eq. (3.14), in powers of the order parameter. Since gauge invariance is not broken in the reference 
ensemble, only even powers of 5" appear. The coefficients in this expansion, i.e. the vertex functions of the effective 
field theory, are connected correlation functions of the anomalous or Cooper channel density ric, Eq. (3.4c), in the 
reference ensemble. In particular, the Gaussian vertex is determined by the pair propagator or anomalous density- 
density correlation fimction in the reference ensemble. Denoting the latter by C{q), the Gaussian term in the LGW 
functional reads 



= E/'^'? [Vir^^^l - C{q)] *"(g) 



(5.1a) 



where we have scaled ^ with l/F^^'. C{q) is the connected propagator C{q) ~ {{ndq) nc{q))o}dis, where {. . .}dis 
denotes the disorder average. C{q) is a complicated correlation function. However, since the reference ensemble is 
by construction a Fermi liquid, the structure of this correlation function is known. Renormalization group arguments 
show that the structure of C at low frequencies and long wavelengths in the limit T ^ is [p2| 



in{n„/{D<f + \n,,\)) 



(5.1b) 



Here flo = ksTo /h is a frequency cutoff on the order of the Debye frequency (for phonon-mediated superconductivity) , 
and Skc is the repulsive interaction in the Cooper channel that is generated within perturbation theory even though 
the bare r^*^-* vanishes. D is the diffusion coefficient of the electrons, h is the renormalized frequency renormalization 
constant whose bare value is H = ttNp/8, see Eq. (3.39a), and Z is the wavefunction renormalization. All of these 
parameters char acter ize the re feren ce ensemble, and it is known that they provide a complete characterization of C{q) 
^Tj. Using Eq. (5.1b) in Eq. (5.1a) we obtain, after again rescaling 



j dq *"*(g) 



1 



ln(l]o/(^q2 + \^n\)) 



(5.1c) 



where t = - {Z\Kc\ - Sk^) /h, with = 7rA^|,r(=)/4. 

Let us now discuss this result for the Gaussian LGW theory. From the structure of the Gaussian vertex we 
immediately read off the values of the exponents r j, 7, a nd z, defined as r(^^(q, O = 0) ~ |qP~'', r(^^(q, rj = 0) ~ , 
and £,T ^ , with the relaxation time, see Eq. (2.1b). They are 



r, = 2 , 7=1 , z = 2 , 
By scaling |q| with the correlation length ^, we also obtain the behavior of the latter. 



(5.2) 



.l/2|t| 



(5.3) 
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The exponent v therefore does not exist, v = oo. 

In order to determine the behavior of the order parameter and the free energy, we need to consider the higher order 
terms in the Landau expansion. The coefficient of the quartic term is a nonhnear anomalous density susceptibihty in 
the reference ensemble which we denote by C^^-* . Due to the cumulant expansion with respect to the disorder average 
there are two different contributions to this coefficient, C^^-* = c|^^ + C'^\ where c|^^ is the disorder average of the 
four-point correlatiau function for a given disorder realization, while C'^'^ is the disorder average of the two-point 
function squared, cj Given that the constant contribution to the Gaussian coefficient, C{q), is barely the leading 
term in the limit q —> 0, and that the quartic coefficient in the magnetic case was singular (see the discussion after 
Eq. (^^)), one docs not expect C^^^ to be finite in the limit of vanishing frequencies and wavenumbers. Indeed, a 
calculation shows that both C^^'' and C2^'' are singular in this limit. Cutting off the singularity by means of a wave 
number |p|, one finds for th^Jeading contributions C'['^'' ^ U4/|p|'* In^ |p| and Cj^-* ^ f4/|p|^~'', respectively, with U4 
and W4 finite coefficicents. cB The same method shows that the most divergent contribution to the coefficient of the 
term of order j^'p" diverges like 

C(2») "f" , (5.4) 

|p|4(n-l) ln2„ ' y ' 

with U2n a finite coefficient. This implies that the Landau expansion of the cutoff regularized LGW theory is an 
expansion in powers of ^f/p^ ln(l/|p|). The order parameter field theory, rather than having a simple LGW form, is 
thus again strongly non-local. This is qualitatively the same effect as in the magnetic case, but here the singularities 
one encounters are even stronger. 

As in the magnetic case, the functional $ can now be analyzed by using standard techniques p^]. We are looking 



for a fixed point where the functional dependence of the 2-point vertex on q and f2„, Eq. (5.1c), is not renormalized. 
This fixes the exponents 77 and z. As in Sees. || and we define the scale dimension of the correlation length to 
be [£] = —1. Power counting then shows that the coefficients U2n {n > 2) of the non-Gaussian terms have scale 
dimensions [u2n] ~ {n— 1)(2 — d), and hence are irrelevant operators with respect to the fixed point for all dimensions 
d > 2. [V4] — —2{d — 2) for 2 < d < 4, and the higher cumulants are even more irrelevant. The upper critical 
dimension is therefore d+ = 2, and for d > 2 the critical behavior obtained from the Gaussian theory is exact. 

The scaling behavior of the order parameter, and of the free energy, is determined by the term of 0{'^^) which is 
a dangerous irrelevant variable with respect to these quantities. For scaling purposes, the cutoff wavenumber |p| can 
be replaced by the inverse correlation length, |p| ~ The scaling behavior "if ~ p^ln(l/|p|) observed above then 
immediately leads to 

, (5.5) 

r I 

One must notice, however, that the order parameter field \1/ is distinct from the physical gap function A. The latter 
determines the gap in the single-particle excitation spectrum, and hence scales like the frequency or like 

A- e-i/l*l . (5.6) 



From Eqs. (5.3) and (5.6) we see that even though the exponents z/ and (3 do not exist, we can assign a value to their 
ratio, 

(3/iy = 2 , (5.7) 

in the sense that A^^ ~ const. 

We next consider the critical behavior of the penetration depth A and the upper critical field ■ Since we 
have shown that the mean-field/ Gaussian theory yields the exact critical behavior at T = 0, all relations between 
observables that are derived within BCS theory are valid. In particular, we have A ^ 1/\/A, and ^ A [5^ ], which 
in conjunction with A ~ yields 



■^''The same structure apppears in the ferromagnetic LGW functional. Here, however, Cj*' is more important than its analog 
in the magnetic case, as we will see below. 

^'^This result for C^/^ holds for 2 < d < 4, and we have neglected possible logarithmic corrections to the power law. For d > 4, 
Cj*' is finite. 
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A'-e , ffc. • (5.8) 

Similarly, we can determine the scale dimension of the conductivity a or the resistivity p. The real part of the 
frequency dependent conductivity in the superconducting phase has a singular contribution, which within BCS theory 
is given by 

Rea,(r!) = y (T„ A,5(r!) , (5.9) 

with On the conductivity in the normal state. For scaling purposes, (5(f2) l/fi ~ 1/A. (t„ is determined entirely by 
properties of the reference ensemble, and hence it does not show any critical behavior and its scale dimension is zero. 
We conclude that the scale dimension of Gs vanishes. If we assume that the conductivity has only one scaling part, 
then the same is true for the conductivity or resistivity in general, and we obtain 

p(t,r) = p(tln5,T62) . (5.10) 

The resistivity thus shows a step discontinuity from a finite value to zero as one crosses the phase boundary at T = 0. 

We now turn to the free energy density /. Hyperscaling suggests that / scales like / ~ T jV ^ ^-('^+2)^ which leads 
to a homogeneity law 

f{t,T,U4) = b~'^'^+^\f{tlTib,Tb^,U4b^-'^) , (5.11) 

where of the irrelevant operators we have written only U4 explicitly. / is proportional to U4 A** ^ l/w4, and hence the 
effective scale dimension of / is [/] = 4, which yields 

f{t,T) = b-^f{t\nb,Tb^) , (5.12) 

Hyperscaling is violated by the usual mechanism that is operative above an upper critical dimension, viz. by means 
of the quartic coefficient being dangerously irrelevant with respect to the free energy. The more exotic violation 
mechanisms that we encountered in the magnetic case are not realized here. By differentiating twice with respect to 
T one obtains the specific heat coefficient ^{t,T), 

j{t,T) = j{t\iib,Tb'^) . (5.13) 

This implies, among other things, that the specific heat coefficient approaches a constant as the temperature is lowered 
to zero at the critical coupling strength. More detailed calculations show a step discontinuity in 7(T = 0) at t = 0. 

The quartic term whose coefficient is 6*2^'' yields corrections to scaling. represents fluctuations in the position of 
the critical point: By making the coupling constant F'"^) a random variable, and integrating out that randomness, one 
obtains a term of that structure. By repeating the arguments of Ref. p9| , one finds that the relative fluctuations of the 
position of the critical point decay anomalously slowly, viz. At/t ~ ^^^-2) ^ Translating that into the corresponding 
fluctuations of Tc via Tc — Tq exp(— 1/|<|) we obtain 

ATjT,^T^^'-^')/^F{\n{To/T,)) , (5.14) 

with F(x) a function that depends parametrically on the disorder. Since it depends on only logarithmically, the 
latter dependence is weak. Equation (5.14) predicts very strong disorder fluctuations in thin superconducting fllms. 
This prediction is in semi-quantitative agreement with recent experiments by Hsu et al. pffl. 



VI. METAL-INSULATOR TRANSITIONS 



A. Disordered Fermi liquids 

In recent years there has been a considerable amount of work on using renormalization group ideas to derive and 
justify Landau's Fermi liquid theory ]6l]] . The main motivation was the idea that non-Fermi liquid notions might 
be important for understanding, for example, the normal state of the high Tc materials, and that in order to derive 
and discuss non-Fermi liquid behavior, one must first have a deeper understanding of why most systems are Fermi 
liquids. In this context, renormalization group ideas are used to understand an entire phase and not just a single 
(although important) critical point in the phase diagram. Technically this implies that one must look for completely 
stable renormalization group fixed points, as opposed to the more usual critical fixed points that are associated with 
phase transitions. 

In this subsection we briefly review how renormalization group ideas can be applied to derive and justify the 
existence of a disordered Fermi liquid phase for d > 2, and, in addition, can be used to derive what are commonly 
known as weak localization effects [|5|. 



29 



1. The disordered Fermi liquid fixed point 



To proceed we parametrize the Q in Eq. ( |3.39b ) in analogy with the parameterization of given by Eqs. ( 3.23 ) 
We write Q in a block matrix form as, 



_ ( yjl - qq'! q 



(6.1) 



where the matrix q has elements whose frequency labels are restricted to n > 0,m < 0. Now, while the sigma- 



model, Eq. (3.39a), can be expressed entirely in terms of q, the corrections to the signia-model action also depend 
on 

AP = P- (P) . (6.2) 

Considering the total action, 

A = ANLaM + AA , (6.3) 

we perform a momentum shell renormalization group procedure. For the rescaling part of this transformation, we 
need to assign scale dimensions to the soft field q, and to the massive field AP as well. If the scale dimension of a 
length L is [L] = —1, we write, in analogy to scaling near a critical point, 

b(x)] = i(d-2 + ry') , (6.4a) 

[AP(x)] = i(d-2 + 77) , (6.4b) 
which defines the exponents rj and rj' . The stable Fermi-liquid fixed point of the theory is characterized by the choice 

r] = 2 , 7?' = . (6.5a) 

Physically, rj' — corresponds to diffusive correlations of q, and rj = 2 means that the correlations of the AP are of 
short range. This is indeed what one expects in a disordered Fermi liquid. In addition, we must specify the scale 
dimension of frequency and temperature, i.e., the dynamical scaling exponent z — [u>] = [T]. In order for the fixed 
point to be consistent with diffusion, that is with frequencies that scale like the squares of wavenumbers, we must 
choose 

z = 2 . (6.5b) 



Now we expand the sigma-model action, Eq. ( 3.39a ), in powers of q. In a symbolic notation that leaves out 



everything not needed for power counting purposes, we write, 

ANLaM = ^ J dxiVqf+H J djcujq^+TT J dxq^ + 0{VWujq\Tq^) , (6.6) 

with the bare coupling constants G ^ I/ctq and H ~ Np. T can stand for any of the three interaction constants r*^^\ 
r^*\ or r'^'^^ Power counting gives 

[G]^[H]^[T]^0 . (6.7) 

These terms therefore make up part of the fixed point action. 

Next consider the corrections that arise within the sigma-model. The leading ones are indicated in Eq. (|6.6|). 
We denote the corresponding coupling constants by Cv2q4, etc., with a subscript that identifies the structure of the 
respective contribution to the action. One finds, 

[C^^,A = [C^,4^[iCT,^f]^~id-2) . (6.8) 



In the last equality in Eq. (6.8) we have considered the square of Cxqi , since any contributions to physical correlation 



functions contain the corresponding term squared. We see that all of these operators are irrelevant with respect to 
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the disordered Fermi-Hquid fixed point for all d > 2, and that they become marginal in d = 2 and relevant in d < 2. 
All other terms in the sigma-model are at least as irrelevant as those considered above. 



It is easy to show that in Eq. (6.3) also contributes to the disordered Fermi-liquid fixed point action. We 
denote this contribution by Ay^*. An inspection shows that AA* depends only on AP and does not couple to q. 
Further, all of the corrections to A ^ ar e more irrelevant for 2 < rf < 4 than the operators with scale dimensions equal 



to 2 — d that are indicated in Eq. ( |6.8D . The scale dimension of the least irrelevant additional corrections is equal to 
-2. 

The above renormahzation group arguments show that the theory contains a disordered Fermi-liquid fixed point 
that is stable for all d > 2. The effective fixed point action is, 

A:ff=ANLaM[q]+AA*[AP] , (6.9) 

and the leading irrelevant terms near this fixed point have scale dimensions given by Eq. (IfL 



2. Scaling behavior of observables 

We now discuss the physical meaning of the corrections to scaling induced by the irrelevant operators that we have 
identified above. Let us denote by the generic name u any of the least irrelevant operators whose scale dimension is 
[u] — — (d — 2), and let us discuss various observables, viz. the conductivity a, the specific heat coefficient jv, the 
single-particle density of states N, and the spin susceptibility Xs- Which of the various operators with scale dimension 
— (d — 2) is the important one depends on the quantity under consideration. 

Le t us first consider the dynamical conductivity, cr(w). Its bare value is proportional to l/G, and according to Eq. 
( |6.7[ ) its scale dimension is zero. We therefore have the scaling law, 

a{uj,u) ^ a{Lub'',ub-^'^-^'>) , (6.10a) 

where b is an arbitrary renormahzation group scale factor. By putting b = 1/w^/^, and using z = 2, Eq. (|6.5l| ), 
as well as the fact that (j{l,x) is an analytic function of x, we find that the conductivity has a singularity at zero 
frequency, or a long-time tail, of the form 

a{Lj) - const + . (6.10b) 

This nonanalyticity is well known from perturbation theory for both noninteracting and interacting electrons |^ , ^ . 
The above analysis proves that the the exact leading nonanalytic behavior. 



The specific heat coefficient, — cy/T, is proportional to the quasiparticle density of states 7? 31 6^, whose 
scale dimension vanishes according to Eq. ( Jf).T\ ) . We thus have a scaling law 

7^(T,u) = 7y(W,u6-(''-2)) , (6.11a) 

which leads to a low-temperature behavior 

7y(T) - const + r('i-2)/2 (6.11b) 

From perturbation theory it is known that shows this behavior only for interacting electrons, while for 
non-interacting systems the prefactor of the nonanalyticity vanishes. This can not be seen by our simple counting 
arguments. 

The single-particle density of states, N , is proportional to the expectation value of Q, and to study the leading 
correction to the finite fixed point valueof N it suffices to replace Q by Q. Then we have, in symbolic notation, 
TV - 1 + (ggt) + . . . = 1 + A7V. The scale dimension of AA^ is [AA^] = 2[g] = d — 2. We find the scaling law 

AN{uj) ^b-^"^-^^ AN{ujb') , (6.12a) 

which leads to the so-called Coulomb anomaly |63| , 

Af(w) - const -I- ^('^-2^/2 . (6.12b) 

Again, this behavior is known to occur only in the presence of electron-electron interactions. 

Finally, we consider the static, wave vector dependent spin susceptibility, Xo(q)- Xo is given by a Q-Q correlation 
function, and the leading correction to the finite Fermi-liquid value is obtained by replacing both of the Q by q. Then 
we have a term of the structure xo ^ T J dx q^q, with scale dimension [xo] = 0. The relevant scaling law is 
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Xo(q,«)=Xo(q,"&-<^-'') , (6.13a) 
which leads to a nonanalytic dependence on the wave number, 

Xo(q) ^const- |q|(^-2) . (6.13b) 

This b ehavior is also known from perturbation theory [ pT[ , and holds only for interacting electrons. As we have seen 



in Sec. IV A 1 above, it has interesting consequences for the theory of ferromagnetism. 

To summarize, we see from the above arguments that all of the so-called weak-localization effects, i.e. nonanalytic 
dependencies of various observables on frequency, temperature, or wave number, in disordered electron systems that 
are well known from perturbation theory, emerge naturally in the present context as the leading corrections to scaling 
near the Fermi-liquid fixed point of a general field theory for disordered interacting electrons. Apart from providing 
an aesthetic, unifying, and very simple explanation for these effects, our arguments also prove that they do indeed 
constitute the leading nonanalytic behavior, a conclusion that cannot be drawn from perturbation theory alone. 

We finally note that nonanalyticities that are very similar to those discussed above occur in classical fluids. In 
that context they are known as long-time tail effects, and they were first discussed theoretically by using many-body 
perturbation theory and mode coupling theory Later, they were examined using renormalization group ideas, 

and they were shown to be related to corrections to the scaling behavior near a hydrodynamic fixed point | |65| . 

B. The Anderson— Mott transition 

It is well known that at sufficiently large disorder, the metallic disordered Fermi liquid phase discussed in the last 
subsection becomes unstable against an insulating phase. Such metal-insulator transitions (MITs) are observed in 
dope d sem iconductors and other disordered electron systems, and the generalized nonlinear sigma-model shown in 
Eq. ( 3.39a| ) is capable of describing them |3^. Metal-insulator transitions whose critical behavior is determined by 



both the electron-electron interaction and the disorder, are commonly referred to as Anderson-Mott transitions, to 
distinguish them from purely disorder driven MITs ('Anderson transitions') and purely correlation driven ones ('Mott 
transitions'), respectively. They can be grouped into two broad classes: (1) Those that are related to fixed points in 
the vicinity oi d — 2, and (2) one related to a Gaussian fixed point in high {d > 6) dimensions. Somewhat strangely, 
these transitions seem to be quite different in nature, and attempts to extrapolate down from d = 6 to d = 3 yields 
results that are incompatible with attempts to extapolate up from d = 2 to d = 3. As a result, the theoretical 
description of the MITs that are observed in 3-d systems is actually an open problem, the substantial amount of 
effort that has gone into this subject notwithstanding. We therefore do not pretend here to be able to give a coherent 
theoretical picture. Rather, we present the results near d — 2 and the results in high dimensions separately, and then 
add some speculations about the behavior in d = 3, with emphasis on the points that are not understood. 

1. Anderson-Mott transition near two dimensions 

This subject has been reviewed in Ref. |Q, and we refer the reader there for details and references to the original 
work. Here we give a brief, updated overview designed to tie this subject in with the more recent developments that 
are reviewed elsewhere in this chapter. 

The MITs near d ^ 2 fall into several distinct universality classes that are related to external fields that couple to 



various degrees of freedom in the spin quaternion space defined in Sec. Ill C 2. Specifically, an external magnetic field 



(to be referred to as MF), magnetic impurities (MI), and spin-orbit scattering (SO) all break some of the additional 



symmetries mentioned after Eq. (3.33) that lead to Q-matrix sec tors I.Q with j,r ^ being soft modes The 
r = 0, z = sector is directly controlled by the Ward identity, Eq. (3.33), and it always is soft. The universality class 



with no additional symmetry breakers is called the generic class, and is denoted by G. Table ^ lists the remaining soft 
modes for all of these universality classes. 

In addition, the critical properties turn out to depend on the type of electron-electron interaction considered; they 
lOije different depending on whether a short-ranged model interaction is used, or a long-range Coulomb interaction. 
c3 For each of the four universality classes G, MF, MI, and SO one therefore distinguishes between short-ranged (SR) 



^^The short-ranged case is not entirely academic. It has been proposed to realize it by putting a grounded metallic plate 
behind the sample which cuts off the long-ranged Coulomb interaction by means of the image charges induced in the plate |68| . 
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TABLE I. Universality classes and soft modes for the MIT near d = 2. 



Universality class Symmetry breaker Soft modes 



MF magnetic field r = 0, 3 ; i = 0, 3 

MI magnetic impurities r = 0, 3 ; i = 

SO spin-orbit scattering r = 0, 1, 2, 3 ; i = 

G none all 



TABLE II. Values for the three independent exponents v, /3, and z for the eight universality classes in d = 2 + e. Values are 
given for d = 2 + e dimensions except for class G, where approximate values for d = 3 based on a two-loop approximation in 
Ref. ]6q| are shown, and for class SO (SR), where the value shown for i/ is the numerical result given in Ref. 0. /3 for the 
class MF (SR) is known to one-loop order, but depends on non-universal quantities, see Ref. The class G (SR) has never 
been considered. 
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z 
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d 


d 


2 + 0(e2) 


? 


fa 5.91 



and long-ranged (LR) subclasses, leading to a total of eight distinct universality classes. For the three independent 
critical exponents (cf. Sec. ||) we choose the correlation length exponent ly, the dynamical exponent z, and the 
exponent P that describes the vanishing of the density of states at the Fermi surface, N^ep) ~ ■ Their values for 
the various universality classes are shown in Table ||. Other critical exponents of interest are related to these three 
by means of exponent relations that are listed in Table H. Which quantities are critical and which are not depends 



on which modes are soft, and therefore on the universality class. This is also indicated in Table III 



The theoretical analysis of t hese universality classes near d = 2 starts from the nonlinear sigma-model, Eq. ( 3.39a ) , 
and the parametrization, Eq. ( |6.l| ) , of the Q~ma'trix or a n equiv alent parametrization. This procedure does not make 



use of the order parameter field theory derived in Sec. Ill C 4, and the fact that the density of states is the order 



parameter for the transition does not become apparent. Rather, one constructs a loop expansion that is equivalent 
to an expansion in powers of q, and proceeds with a rcnormalization group analysis analogous to that for the 0{N)- 
Heisenberg model in d = 2 + e dimension |2^. For the universality classes MI, MF, and SO, this procedure readily 
leads to fixed points corresponding to MITs, and the asymptotic critical properties at these transitions have been 
worked out to lowest order in a 2 + e expansion p5[ . For the conductivity and the density of states, the results are in 
reasonable agreement with experiments, although the ill-behaved nature of the 2 + e expansion precludes quantitative 
comparisons Measurements of the spin susceptibility and the specific heat, on the other hand, have shown 

non-Fermi liquid like behavior on either side of the MIT [p9[ , which the nonlinear sigma-model can not explain. Very 
similar results were obtained experimentally on systems that are believed to be in the universality class G |^0| . These 
features are usually interpreted in terms of local magnetic moments that are believed to be ubiquitous in disordered 
electronic systems |7l| . Currently no unified theory exists that can describe the interplay between these local moments 
and the critical behavior of the transport properties at an MIT. 

For the universality classes MI (SR) and SO (SR) it turns out that the electron-electron interaction amplitude scales 
to zero under renormalization p6| . The critical behavior in these cases is therefore the same as for the corresponding 
universality classes of non-interacting electrons |43H7^ . The density of states is uncritical (/3 = 0) , and the dynamical 
exponent is z = d. For the class MI (SR) one finds a fixed point at one-loop order that corresponds to a MIT, with 
a correlation length exponent v as given in Table |l|. The behavior for the class SO (SR) was unclear for a long time, 
since in this case there is no fixed point that is accessible by means of a 2 + e expansion. However, there is now good 
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TABLE III. Physical quantities, their critical behavior, and relevant exponent relations. Also listed are the universality 
classes for which the respective quantities show critical behavior. The order parameter density is the density of states. The 
first two columns are generally valid, while the last two pertain to the Anderson-Mott tran sition near d — 2 only. In high 
dimensions, there is only one universality class, and some exponent relations change, see Sec. VI B 3. 



Physical quantity 



Scaling behavior 



Exponent relation 



Universality classes 



Correlation length 

Correlation time 

Density of states 

Electrical conductivity 

Specific heat 

Heat diffusion coefficient 

Spin diffusion coefficient 

Order parameter density 
correlation function 

Order parameter density 
susceptibility 



N(€F,t) ~ 

a{t,T = = 0) ~ 
cv{t = 0,T) ^T^-'"/^ 
Dh{t,T = 0) f"^ 
Ds{t,T = 0) 
G(k,t = T = 0) - 

Xop{t,T = 0)^t-^ 



independent exponent 
independent exponent 
independent exponent 

s = v{d~2) 

K — z — d 
Sh = i^{z - 2) 
Ss — v(z — 2) 
r) = 2 -'d + 2f3/u 

7 = i/d - 2/3 



all 
all 

MI(LR), MP, SO(LR), G 
all 

MI(LR), SO(LR) 
MI(LR), SO(LR), G 
G 

MI(LR), MP, SO(LR), G 
MI(LR), MP, SO(LR), G 



evidence from both numerical analyses and high-order perturbation theory for a MIT in c? = 3 |^,^,^ . 

For the universality classes MI (LR), SO (LR), and MF, the renormalization group analysis yields MIT fixed 
points where the interaction strength scales to a finite number, so that the critical behavior is different from that of 
noninteracting electrons. The results are listed in Table ||. For the MI and MF universality classes the Cooper channel 
(r = 1,2 in the spin-quaternion basis) does not contribute to the soft modes, see Table |. It is therefore sufficient 
to rcnormalize the particle-hole degrees of freedom. For the class SO (LR), on the other hand, one has to deal with 
the Cooper channel as well. How to properly do this in an interacting system has given rise to some controversy. 
The current state of affairs is that no satisfactory renormalization group treatment of the Cooper channel interaction 
amplitude F^'^^ exists |Q. However, there is general agreement that F^'^' is irrelevant at the MIT fixed point, so that 
the asymptotic critical behavior is unaffected by this unsolved problem. It is also likely that F'-'^^ is only marginally 
irrelevant, and gives rise to logarithmic corrections to scaling ||7q , |75|| . 

The behavior one encounters in class G is substantially more complicated than that in any of the other universality 
classes. Low-order applications of the perturbative renormalization group reveal a runaway behavior of the spin- 
triplet interaction constant, F^*^ and the flow equations do not allow for a fixed point corresponding to a MIT [ p4| . 
There is a popular belief that this runaway flow is an indication of local-moment formation, which weak-disorder 
perturbation theory is inadequate to describe [Q. We do not believe that this is the correct interpretation of the 
runaway flow, for the following reasons. As has been shown in Ref. ||7^, the leading singularities in the runaway 
region, i.e. for the renormalized F^'^ — > oo, can be resummed exactly to all orders in the loop expansion. The result 
is a pair of coupled integral equations for the spin and heat diffusion coefficients, the solution of which yields a phase 
transition where the spin diffusion coefficient vanishes. This transition has also been studied with renormalization 
group techniques, which corroborated the results obtained from the integral equations | |78[ . The conductivity or charge 
diffusion coefficient was found to be uncritical, and decoupled from the spin transport in the critical region, so the 
transition is not a MIT. A very remarkable aspect of these results was that the critical behavior could be obtained 
exactly, apart from possible logarithmic corrections, in all dimensions d > 2. The reasons for this unusual feature 
were not clear at the time. Originally, the nature of this transition was not known, and the original papers speculated 
about it leading to a pseudomagnetic phase with 'incompletely frozen spins'. While such an 'incompletely frozen spin 
phase' might have had something to do with local moments, later developments showed that these speculations were 
incorrect. The exact critical behavior found in Refs. ||7^,Q turns out to be identical with that of the exactly soluble 
paramagnet-to-ferromagnet transition in disordered itinerant electron systems that was reviewed in Sec. IV A 1 above. 
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G 

FIG. 4. Schematic phase diagram for disordered itinerant electron systems at T = close to d = 2, in the plane spanned 
by the spin-triplet interaction constant F*-'' and the disorder G. The phases shown are the paramagnetic metal (PM), the 
ferromagnetic metal (FM), and the insulator (I). It is not known whether there is another phase transition within I from a 
ferromagnetic to a paramagnetic insulator. 



This can hardly be a coincidence, and we therefore beheve that the correct interpretation of the runaway flow found 
in the G universaUty class is a transition to a long-range ordered ferromagnetic phase, not an instability against 
local moment formation. The fact that the runaway flow occurs for arbitrarily small values of the disorder and the 
interaction strength, provided that one works close to d = 2, suggests a phase diagram for low-dimensional systems 
as shown in Fig. The point is that close to two dimensions, the paramagnetic metal phase and the insulator phase 
are always separated by a ferromagnetic phase. This is the reason why the MIT in the G universality class is not 
accessible by means of an e-expansion about d — 2, and it is what is causing the runaway flow in the 2 + e expansion. 
Whether or not the paramagnetic metal phase is unstable against the formation of local moments, and if so, why this 
physical feature is not reflected in the sigma-model, is a separate question. We also point out that there are other 
competing instabilities in this universality class, in addition to the ones noted above. For example, in Ref. ||79|] we 
have argued that in d = 2, there is an instability to a novel type of even-parity superconductivity. All of this suggest 
that the actual phase diagram in and near two dimensions is very complicated and, at this point, not understood. 

For d > 2 + Cc, with ec a critical value that was estimated in Ref. ||7^, there is a direct transition from the 
paramagnetic metal to an insulating phase, see Fig. ^. Obviously, this transition cannot be described by means of a 
controlled e-expansion about d = 2. The properties of this transition in d = 3 have been estimated by means of a 
two-loop approximation , the results of which are shown in Table ||. 

Finally, we mention that all existing treatments of the Anderson-Mott transition conclude that d~ = 2 is the lower 
critical dimension of the problem, in the sense that in d = 2 there is no metallic phase for any degree of disorder. Until 
recently, the experimental situation was in agreement with this. Experiments that seem to show a transition from 
an insulating phase to a conducting one in 2-d Si MOSFET systems therefore came as a considerable surprise [|oj. 
Several theoretical speculations have been put forward |]8l|l , but the nature of this phenomenon is currently unclear. 



2. Anderson-Mott transition in high dimensions 



Let us now return to the generalized sigma-model in the form given by Eq. (3.42), where the soft modes have 
been integrated out. This rewriting of the action one would assume to lead to a simple formulation of the problem 
in sufflciently high dimensions, where the soft modes do not provide the leading physical effects. This expectation 
will prove to be correct, and will allow us to determine the exact critical behavior at a metal-insulator transition in 
dimensions d > 6 . 

We look for saddle-point solutions Qsp, A^p that are spatially uniform and satisfy 
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G 

FIG. 5. Schematic phase diagram for a 3-d disordered itinerant electron system in the plane spanned by the spin-triplet 
interaction constant F'*' and the disorder G. The phases shown are a paramagnetic metal (PM), a ferromagnetic metal (FM), 
a paramagnetic insulator (PI), and a ferromagnetic insulator (FI). M denotes a multicritical point. 



.14a) 



6 ai'^^^ 



.14b) 



where the superscript (0) denotes the saddle-point approximation. The replica, frequency, and spin-quaternion 
structures in Eqs. (3.14) are motivated by the fact that (^Q^m) ^'^'^ (r^nm) have these properties, and that in the 
mean-field approximation averages are replaced by the corresponding saddle-point values. 
Taking the extremum of the action, Eq. ( ^.42 ), with respect to A or one obtains 



(^r)' = i-Gr(^)/„(^(°)) 

and taking the extremum with respect to Q or N'^^^ gives 



.15a) 



.15b) 



/„ in Eq. ( |6.15a ) is a functional of ^ that contains several frequency integrations over S^} , and that in addition depends 



parametrically on the external frequency w„, on the disorder, and on the interaction constants | |38| . For simplicity 
we show only the spin-singlet interaction, F^^^, explicitly. It has been shown in Ref. |Q that the other interaction 
channels do not qualitatively change the conclusions. 



Nn is simply related to the density of states through Eq. ( 3.31 ). For fixed F^'*), Eqs. ( 6.15 ) describe a decrease of 
the density of states with increasing disorder G. Furthermore, N{eF +oj) is nonanalytic aX lj = Q. An it eration to first 
order in G recovers the 'Coulomb anomaly' that is well known from perturbation theory |^3|, see Eq. (6. 12b) above. 



With further increasing disorder, N{ep) vanishes at a critical value Gc of the disorder. The mean-field approximation 
thus describes a phase transition with a vanishing density of states at the Fermi level, the hallmark of the Anderson- 
Mott transition. (Transport properties we will discuss shortly). Note that in the absence of interactions, iV^^^ = 1. 
This refiects the fact that the density of states does not vanish at an Anderson transition . 

The saddle point equations, Eqs. (|6.15|) , together with an expansion to Gaussian order about the saddle point, 
determine the critical behavior of this Anderson-Mott transition in high dimensions. One finds standard mean- 
field/Gaussian values for all static exponents. 



7 = 1 



r; = 



(5 = 3 



(6.16a) 
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and for the dynamical critical exponent, 



z = 3 . (6.16b) 
The mean-field critical behaviors of the charge and spin diffusion coefficients have also been determined Bq], with 



the result that they both are proportional to the order parameter. For the exponents s and Sg defined in Table III 

.s E 



this means 

s = s, = 1/2 (6.16c) 

We now must ask what the upper critical dimensionality for the problem is, i.e. the dimension djj" above which the 
above critical behavior is exact. A superficial inspection of the action suggests = 4. However, it turns out that 
under renormalization additional contributions are created that were not in the bare action, and that lead to (i+ = 6 
psf . The structure of these additional terms is reminiscient of the structure encountered in the theory of magnets in 
random magnetic fields. This becomes plausible if one recalls that the random potential in the action couples to the 



electron density, see Eq. (3.2c). If we Fourier transform from imaginary time to Matsubara frequencies, we obtain a 
term proportional to 

/ (ixu(x)^?/5„(x)'(/'„(x) . (6.17) 

Since the expectation value (■0„(x)'(/'n(x)) determines the density of states, this means that the random potential 
couples to the order parameter field for the phase transition, just like a random magnetic field couples to the magnetic 
order parameter. A detailed technical analysis does indeed confirm this analogy. Not only is c?+ = 6, but many 
structures of the theory coincide with those encountered in the theory of the random-field Ising model , including 
a quartic coupling constant, u, that acts as a dangerous irrelevant variable even below the upper critical dimension 
= 6. An explicit 6 — e expansion of the critical exponents yields also results that, at least to one-loop order, 
coincide with those for the random-field Ising model. 

The implications of these results for the critical behavior in d = 3 are discussed in the next subsection. For 
dimensions d > 6, the mean-field results reviewed above represent the first example of an Anderson-Mott transition, 
or metal-insulator transition of disordered interacting electrons, for which the critical behavior has been determined 
exactly. Remarkably, this transition has nothing obvious in common with the Anderson-Mott transitions that are 



described by the same model in rf = 2 + e and that have been discussed in Sec. VI B 1 above. In particular, there is no 
sign here of the various universality classes that are present in d = 2 + e. This is not just the strong universality that 
one expects from a mean-field theory, as it persists in the 6 — e expansion. It is therefore likely that the Anderson-Mott 
transition described above is of a different nature than those found in d = 2 + e. Which description is closer to what 
actually happens in 3-d systems is currently not known. 



3. Anderson-Mott transitions in three dimensions: Conventional scaling scenario 



The similarities between the Anderson-Mott transition in high dimensionalities and random-field magnets men- 
tioned in the last subsection gives rise to a scaling description of the former p8|| that is quite different from the one 
presented in Sec. 



VI Bl 



It is important to point out that for the case of random-field magnets, this scaling scenario 
almost certainly is not correct. However, there are subtle technical differences between the two systems that suggest 
that the situation may be different in the case of the Anderson-Mott transition. We therefore present this scenario, 
which is a conventional (i.e., power-law) scaling description, here, and a different, more exotic, possibility in Sec. 
VI B 4 below. Which of these two possibilities provides a better description of the actual behavior in d = 3 is not 



known. 

There are two crucial ingredients to a conventional scaling description of the Anderson-Mott transition based on 
the order parameter formalism of Sees. |III C 4 and VI B 2. The first one is the realization that the electron-electron 
interaction is renormalization group irrelevant at this transition. Efl As a result, the frequency mixing that produced 



■^''The conductivity a is related to the charge diffusion coefficient Dc according to the Einstein relation a = Dcdn/dfj,, with 
dn/dfi the thermodynamic compressibility. The result quoted for s assumes that dn/dn has a noncr itical contribution. 
^*This is true despite the fact that the interaction is necessary for the transition to occur, see Eq. (6.15a). 
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nontrivial dynamical exponents z in the 2 + e expansion, is absent here, the frequency plays the role of an external 
field that is conjugate to the order parameter, and z is not independent. If we denote the scale dimension of the 
external field by yh-, we have 



z^Vh^ (5/3/ V 



(6.18) 



The second important point is the existence of a quartic coupling constant, u, that acts as a dangerous irrelevant 
variable with a scale dimension of [u] = —d. This adds a third independent static exponent, 0, to the usual two. The 
order parameter obeys a scaling or homogeneity relation. 



N(t, n,u,...) = 6(2-<i-'))/2 7v(t&i/^ . . .) 



(6.19a) 



Here we have denoted by z the scale dimension of or T, [fl] = [T] = z. Upon elimination of the dangerous irrelevant 
variable u, z turns into the effective dynamical exponent z, 

7V(t,0) = 6(2+e-rf-'7)/2 _ (6.19b) 

This relates the OP exponent (3 to the three independent exponents v, rj, and 9 through the scaling law, 

/3=^(d-6l-2 + r?) . (6.20a) 

The remaining static exponents are given by the usual scaling laws, with d —t d — 6 due to the violation of hypcrscaling 
by the dangerous irrelevant variable. 



5={d-6 + 2-ri)v/2l3 , 
7 = 1/(2 - rf) . 



(6.20b) 



Now we consider the thermodynamic susceptibilities dn/d^, 7y, and Xs- From general arguments given in Ref. ]3^ ] 
we expect all of them to share the same critical behavior. Denoting their singular parts collectively by Xsing, we have 



Xsi7ig{t,T) = b 



_ u-d+0+z 



Xs^ng{th^'\Tb') 



(6.21) 



This links the static critical behavior of the thermodynamic susceptibilities, characterized by the exponent k, 
Xsing{t) - (see Table pi, to that of the OP, 



K = (3 



(6.22) 



where we have used the scaling laws, Eqs. (6.18), (6.20b). The fact that all of the thermodynamic sysceptibilities scale 
like the order parameter is a consequence of the random-field or static disorder fluctuations being dominant over the 
quantum fluctuations. 

Let us also consider the scaling behavior of the transport coefficients. The charge, spin, and heat diffusion coeffi- 
cients, which we collectively denote by Z?, all obey the same homogeneity relation. 



By the definition of the exponents Ss and Sh (see Table [II), this yields 



(6.23) 



(6.24) 



If the thermodynamic density s uscep tibility, dn/d u, has a noncritical background contribution, then the conductivity 
exponent s is also given by Eq. {^^ . Using Eq. (16 20a| ), we find the following generalization of Wegner's scaling law, 
Eq. (^lb|) and Table |l|. 



■s^-{d-2-0-7j) 



With these results, we can write the homogeneity law for the static electrical conductivity as 



(6.25) 



(6.26) 



where the scaling function is determined by a(t = l,T), i.e. by the conductivity far from the transition. The 
salient point is that, according to this scaling scenario, a/f^ in the critical region is a function of T/f^^, rather than of 
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FIG. 6. Dynamical scaling plot of the conductivity data from Fig. 1 of Ref. |8^. The plot assumes a critical stress 
Sc = 6.71 kbar, and exponent values s = 0.29, vz — 1.82. Only data in the temperature range T < 60 mK have been 
included in the plot, and different symbols denote different stress values, from 5* = 6.59 kbar to S = 8.03 kbar. We have chosen 
To = 100 K, and the relation between t and S — Sc was taken from Ref. Q, viz. t = {S - Sc) 5.4 x 10"^ (kbar)-\ The inset 
shows that the data cease to scale once the temperature region 60 mK < T < 225 mK is taken into account. (From Ref. Q.) 



T and t separately. This can be checked experimentally. In Eq. ( 6.26| ), the exponent s is related to the independent 
exponents v, rj, and 9 by Eq. ( |6.25 ). An important point in this context is again the Harris criterion which requires 
1/ > 2/ d. If Wegner scaling, Eq. ( ^.4b| ), were valid (as it is, e.g., in the transitions in d = 2 + e discussed in Sec. 
VI B 1), then it would follow that j/ > 2/3 in three dimensions. This puts a severe constraint on any interpretation 
of experimental results, and has been a much-discussed issue in the case of the metal-insulator transition that is 
observed in Si:P A thorough discussion of these experiments in the light of Eq. (5.26) has been given in Ref. [ p4| . 
The net conclusion was that the data allow for reasonable dynamical scaling plots, provided that one chooses values 
of s that are substantially smaller than 2/3. An example of such a scaling plot is shown in Fig. ^ This indicated that 
an interpretation of this experiment in terms of conventional power law scaling is possible only within the framework 
of a theory that, like the present one, violates Wegner scaling. This seems to favor the order parameter description 
of the Anderson-Mott transition over the (2 + e)-dimensional approach that was reviewed in Sec. VIB 1 . However, a 
competing experiment on the same material fs^ obtained results that are incompatible with Ref. |85|J86(| , and allow for 
scaling plots of comparable quality with a value of s « 1.3. This experimental controversy has not been resolved [ p8| , 
and the experimental accuracy is not sufficient for a scaling plot to discriminate among them. For a more detailed 
discussion of this point we refer the reader to Ref. [r3 . 



4- Anderson-Mott transitions in three dimensions: Activated scaling scenario 

In the previous subsection we gave a scaling theory for the Anderson-Mott transition, assuming it was a conventional 
continuous phase transition, albeit with dangerous irrelevant variables playing a crucial role. There is, however, 
another possibility. We have seen that the Anderson-Mott transition near d = 6 has structural similarities to the 
transition that occurs in a classical random-field Ising model near d = 6. The latter has been predicted and observed 
to display glass-like features, and activated rather than power-law scaling in d = 3. If the analogy to the Anderson- 
Mott transition still holds in d = 3 (which at this point is merely a subject of speculation), then the remarkable 
conclusion is that the Anderson-Mott transition also has glassy aspects. This possibility has been explored in some 
detail in Ref. |Q . Here we briefly show how to modify the scaling theory for the Anderson-Mott transition in order 
to allow for this possibility. 

The chief assumption in a scaling theory for a phase transition that shows activated scaling is that the critical time 
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scale grows exponentially with ^, 



\n{T/To)^e , (6.27) 

with To a microscopic time scale, and ip a generalized dynamical exponent. Physically, this equation implies an 
exponential growth of the relaxation time as the transition is approached, i.e., behavior that is typical of a glass 
transition. As a result of this extreme critical slowing down, the system's equilibrium behavior near the transition 
becomes inaccessible for all practical purposes, and the systems falls out of equilibrium on realizable experimental time 
scales. Below we will point out some of the experimental consequences of this. Before we do so, let us briefly discuss 
the physical ideas behind our picture of a 'glassy' Anderson-Mott transition. As already noted, the Anderson-Mott 
transition problem contains an intrinsic frustration feature: Electron-electron interactions always suppress the local 
density of states, i.e. the order parameter, while the random potential can cause either local increases or decreases 
in the order parameter. This frustration means, for example, that local insulating clusters exist inside the metallic 
phase. Elimination of these clusters requires a large free energy barrier to be overcome. These barriers are assumed 
to grow like L"^, w ith L some length scale, which near the critical point is given by ^. Via the Arrhenius law, this 



leads to Eq. (6.27) 



To construct a scaling theory for this type of transition, the usual homogeneity laws need to be generalized to allow 
for activated scaling. Since the barriers are expected to be normally distributed, while the relaxation times have a 
much broader distribution, the natural scaling variables are In t/tq or ln{To/T), with Tq a microscopic temperature 
scale such as the Fermi temperature. Considering a variable Q, which we assume to be self-averaging, we therefore 
expect a homogenity law 

Q(t,T) = 6-t«lF,(t.V.,_|_) , (6.28) 
where Fq is a scaling function. For the tunneling density of states this implies, 

N{t, T) ^ b-^'^FN (th^'" , TT^^-rr) ■ (6.29a) 



in(ro/r) 

Eliminating the parameter h by choosing h"^ = \n{Ti^/T) gives 

m.T)= ^ Gnif'^HT./T)] , (6.29b) 

[ln(To/T)] 

with Gn another scaling function. Equation ( |6.29lj ) illustrates a general and important result: At criticality, t = 0, 
the critical singularities will be only logarithmic in ln(To/T). This is the quantum analog of the fact that classical 
static glass transitions are experimentally inacessible. 

Other quantities that are singular at the critical point have been discussed in Ref. [Q . Here we explicitly consider 
the magnetic susceptibility, and the specific heat coefficient. Given that the entropy and the magnetization, as 



thermodynamic quantities, satisfy homogeneity laws like Eq. (3.28 ), it has been argued [ p4[ that neither Xm nor 7y 
satisfy such a law. As a consequence of this, one finds that these quantities are singular even in the metallic phase. 
For example, Xm scales like 

rp—l + constxt""'' 

(j) is an exponent related to the difference in scaling between T and a magnetic field H. The singular behavior in 
the metallic phase means that there is a Griffiths phase, or a region away from the critical point in which certain 
observables become singular at various values of t. 

Finally, it has been argued in Ref. [Q that the conductivity, a, is not self averaging and that instead one should 
consider l„ = {\n.((jQ/ a)} with ctq the bare conductivity. The physically reason for this is that u is related to a 
relaxation time which, as mentioned above, is log-normally distributed, obeys a scaling law 



Ut,T)^b'f'F^ ( tb^'" , , , , 
V ln(7o/r), 

= ln(To/T) a {f^^ HTo/T)) . (6.31) 
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The important point is that the observable conductivity at T = will vanishes exponentially fast as t ^ 0, but that 
it is impossible to observe this behavior because it requires exponentially low temperatutes. The same arguments also 
imply that one should expect large sample-to-sample variations in the measured conductivity. 

There is some experimental support for the glassy picture presented above. First, it has been known for some 
time that there are large sample-to-sample variations in the conductivity measurements in doped semiconductors. 
Already in the early 1980's this was observed in Si:P for t < 10^^ in the mK temperature range Later, similar 

effects were blamed on some type of thermal decoupling The so-called Griffiths like phase in which Xm and 

"fv are singular already in the metallic phase has also been seen in the doped semiconductors. Coventionally this 
behavior has been attributed to local moment formation |^ , and the vanishing of the Kondo temperature in 
disordered systems. Whether, or how, these effects are related to the Griffiths singularities that arise in our glassy 
phase transition theory is not obvious and warrants further work. 
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